Decoherence in weak localization I: Pauli principle in influence functional 
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This is the first in a series of two papers (I and II) , in which we revisit the problem of decoherence 
in weak localization. The basic challenge addressed in our work is to calculate the decoherence 
of electrons interacting with a quantum-mechanical environment, while taking proper account of 
the Pauli principle. First, we review the usual influence functional approach valid for decoherence 
of electrons due to classical noise, showing along the way how the quantitative accuracy can be 
improved by properly averaging over closed (rather than unrestricted) random walks. We then 
use a heuristic approach to show how the Pauli principle may be incorporated into a path-integral 
description of decoherence in weak localization. This is accomplished by introducing an effective 
modification of the quantum noise spectrum, after which the calculation proceeds in analogy to the 
case of classical noise. Using this simple but efficient method, which is consistent with much more 
laborious diagrammatic calculations, we demonstrate how the Pauli principle serves to suppress 
the decohering effects of quantum fluctuations of the environment, and essentially confirm the 
classic result of Altshuler, Aronov and Khmelnitskii for the energy-averaged decoherence rate, which 
vanishes at zero temperature. Going beyond that, we employ our method to calculate explicitly the 
leading quantum corrections to the classical decoherence rates, and to provide a detailed analysis of 
the energy-dependence of the decoherence rate. The basic idea of our approach is general enough to 
be applicable to decoherence of degenerate Fermi systems in contexts other than weak localization 
as well. — Paper II will provide a more rigorous diagrammatic basis for our results, by rederiving 
them from a Bethe-Salpeter equation for the Cooperon. 



I. INTRODUCTION 

The weak localization of electrons by coherent 
backscattering in a disordered conductor, which mani- 
fests itself via a characteristic contribution to the mag- 
netoconductivity, is a unique, particulary robust inter- 
ference effectiiSiiiiSiSii. It is not suppressed by thermal 
averaging, and the temperature dependence of the effect 
arises only due to the destruction of quantum coherence 
by inelastic scattering events, whose likelihood increases 
with rising temperature. The study of decoherence, and 
in particular of the temperature dependence of the deco- 
herence rate 7 ¥ >(T) governing the magnetoconductivity, 
therefore plays a central role in this subject. 

There are two features which make this problem non- 
trivial, related to the influence of low- and high-frequency 
environmental fluctuations on the propagating electron, 
respectively: On the one hand, the environmental fluc- 
tuations at the lowest frequencies do not contribute to 
decoherence, since they are so slow that they resemble 
an elastic impurity potential: for trajectories of duration 
t, environmental frequencies Lb < 1/t do not contribute 
to decoherence. This fact is most easily accounted for 
in an influence functional or path-integral description in 
the time domain, which works well as long as the fluctu- 
ations are classical. On the other hand, environmental 
modes with frequencies much higher than the tempera- 
ture do not contribute either, since an (electron- or hole- 
like) quasi-particle propagating with energy e ~ T rela- 
tive to the Fermi surface does not have enough energy to 
excite them: Pauli blocking forbids the quasi-particle to 



loose an energy Co larger than ~ T to the environment. 
This fact is obvious in the frequency domain, where Pauli 
factors such as f{e) [1 — f(e — ui)] become explicit (/ be- 
ing the Fermi function); hence, a proper treatment of 
high frequencies is most easily achieved in a perturbative 
many-body calculation in the frequency domain, which 
allows for a fully quantum-mechanical treatment of the 
environment. 

Although the essential physics of both the low- and 
high-frequency environmental modes is well understood, 
it is rather difficult to explicitly and accurately treat 
both regimes on an equal footing within a single, uni- 
fied framework. On the one hand, standard influence 
functional approaches usually do not incorporate the 
Pauli principle explicitly (a notable exception being the 
work of Golubcv and Zaikin^SiSiiSiiiiiSiiiM which is, how- 
ever, controversial 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 and whose re- 
sults for 7 V (T) we disagree with). On the other hand, 
diagrammatic approaches in the present context have dif- 
ficulties in accurately dealing with infrared divergencies, 
which are often simply cut off by hand, with the cutoff 
being determined self-consistently (or else the presence 
of an external cutoff is assumed^, as provided by an 
applied magnetic field) . In the present series of two pa- 
pers (I and II2&), we fill in the respective "gaps" in both 
the influence functional approach (paper I) and the dia- 
grammatic approach (paper II), by showing how each can 
be extended to achieve an accurate, explicit treatment 
of both low- and high-frequency modes. The resulting 
two approaches, though thoroughly different in style and 
detail, yield the same result for the Cooperon decoher- 
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cncc rate 7^(T), for which wc evaluate both the lead- 
ing and next-to-leading terms (in an expansion in which 
the dimensionless conductance is the small parameter). 
The leading terms coincide with that found by Altshulcr, 
Aronov and Khmelnitskii^ 5 . (AAK) for decoherence due 
to the thermal part of Nyquist noise (which we shall 
call classical white Nyquist noise below). The next-to- 
lcading terms are checked against and found to be con- 
sistent with results for the magnetoconductivity in large 
magnetic fields of Aleiner, Altshuler and Gerzhensoni^ 
(AAG). 

Paper I is intended for a wide audience and will hope- 
fully be accessible to nonexperts. It presents a path inte- 
gral analysis of decoherence by quantum Nyquist noise, 
achieving not only a natural infrared cutoff, but also in- 
corporating the Pauli principle in a physically transpar- 
ent way, by suitably modifying the interaction propaga- 
tor. In particular, wc offer an elementary but quantita- 
tively accurate explanation for why and how the Fermi 
function enters the decoherence rate. As a byproduct of 
our analysis we show (i) how the accuracy of the path- 
integral approach can be improved by performing tra- 
jectory averages over closed (as opposed to unrestricted) 
random walks^ 7 -, (ii) calculate the leading quantum cor- 
rections to the classical results for the decoherence rate, 
and (iii) are also able to analyse explicitly the energy 
dependence of the decoherence rate. 

We reach our goal by a scries of steps, whose main 
arguments and results are summarized concisely in Sec- 
tion[n]in a type of overview for the benefit of readers not 
interested in the details of the derivations. 

The price paid for our avoidance of a large formal appa- 
ratus in favor of simple, transparent arguments is that pa- 
per I does not entirely stand on its own feet: its discussion 
of Pauli blocking relies in part on heuristic arguments 
and/or influence functional results derived elsewher e - 9 ! 24 . 
In paper Il2&, addressed to experts, we aim to put the 
heuristic arguments of paper I on a solid footing, by 
redcriving the main results for the Cooperon propaga- 
tion in a completely different manner, using purely di- 
agrammatic means. To this end, we use Keldysh per- 
turbation theory to set up a Bethe-Salpeter equation for 
the Cooperon, which includes both self-energy and ver- 
tex contributions to the Cooperon self-energy and whose 
leading terms are free from both ultraviolet and infrared 
divergencies. This equation is then converted to the time 
domain and solved approximately with an exponential 
Ansatz, which, remarkably, turns out to reproduce the 
results of the paper I. 

Our work is built on a foundation laid over many years 
by many different authors. The influence of classical 
(purely thermal) white Nyquist noise was first studied 
in the seminal work of AAK™, where they derived a 
path-integral description and were able to solve exactly 
the quasi 1-dimensional case. Chakravarty and Schmid 
elaborated on this approach in their review 6 -, which also 
includes a detailed discussion of electron-phonon scatter- 
ing. More recently, Voclkcr and KopietaS 8 . provided an 



alternative to path-integration, an "Eikonal" ansatz for 
the time-evolution of the Cooperon, which also includes 
the correct infrared behaviour. 

All of these works, explicitly or implicitly, deal with 
the Pauli principle by using a classical noise spectrum 
that is derived from the physical quantum-mechanical 
spectrum by eliminating the possibility of spontaneous 
emission into the bath (see our discussion in Section 
IV A|) . This prescription was consistent with pcrturba- 
tive diagrammatic calculations (such as the calculation 
of the inelastic electron scattering rate in Ref. [2^) , and 
it was recently reconfirmed by AAG™ via a detailed di- 
agrammatic calculation of the short-time behaviour of 
the Cooperon, to leading order in the interaction. An 
expansion in the quantum corrections to the picture of 
decoherence by purely classical noise, performed by Vav- 
ilov and Ambegaokar— , yielded similar results. 

These recent studies 1 ^ 7 - were motivated by and con- 
tributed to a controversy that arose when Golubev and 
Zaikin (GZ) claimed 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 to have demonstrated 
theoretically that electron-electron interactions intrinsi- 
cally cause the decoherence rate 7^ to saturate at a 
nonzero value at low temperatures, and that this ex- 
plains the saturation that has been observed in some 
experiments™ In these papers, GZ proposed a new, 
exact Feynman- Vernon influence functional for electrons 
under the influence of an environment, which takes 
proper account of the Pauli principle (as confirmed in 
Ref. I24F) . However, the evaluation of this influence 
functional is not straightforward, and the approxima- 
tions which GZ adopted to this end have been heav- 
ily criticized, 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 , in particular those 
pertaining to the terms associated with Pauli blocking. 
Very recently, von Delft has showr™ that if the Pauli 
blocking terms are treated somewhat more carefully to 
include recoil effects, GZ's approach actually does repro- 
duce the celebrated results of AAK for the decoherence 
rate 7 V (T), which does not saturate at low temperatures. 
The analysis of Ref. |^ constitutes a formal counterpart 
to the present paper I, in which we use partly heuristic 
arguments to reach the same conclusions as Ref. in a 
more intuitive manner. 



II. OVERVIEW OF RESULTS 

Before embarking on a detailed calculation of the de- 
coherence rate, we present in this section an overview of 
the main results and arguments contained in the present 
paper, and a short analysis of their various strengths and 
weaknesses. It is hoped that the reader will thereby gain 
a birdseye view of the problems that typically arise in cal- 
culations of 7^, a feeling for what is required to conquer 
them, and a glimpse of the type of results obtained. 

The weak localization contribution to the magnetocon- 
ductivity of a quasi d-dimensional disordered conductor 
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can be written in the forrn^ 

SaJ h = d - dtC(0,t). (1) 

Here v = TOfcp/(27r 2 ?i 2 ) is the 3-dimensional density of 
states per spin at the Fermi surface, r e i the elastic scat- 
tering time, and Ud = 2e 2 VdD is the sample's classical 
Drude conductivity for d = 3, the inverse square resis- 
tance for a d = 2 film of thickness a, or the inverse re- 
sistance per length for a d = 1 wire of cross sectional 
area a 2 , with Vd — a?~ cl v being the effective density of 
states per spin of the corresponding dimensionality, and 
D = VpT i/3 the diffusion constant. C(r,t) denotes the 
Cooperon propagator, in the position-time representa- 
tion, in the presence of interactions and a magnetic field. 
For r — it gives the probability for an electron prop- 
agating along two time-reversed paths to return within 
time t to the starting point without losing phase coher- 
ence, thus enhancing the backscattering probability and 
reducing the conductivity. In the absence of decoherence 
and a magnetic field it is given by the classical diffusion 
probability density (the "diffuson" ) . 

For ease of reference, our notational conventions will 
mostly follow those used in Ref. 0- I 11 particular, var- 
ious incarnations of the Cooperon propagator will occur 
below, related by Fourier transformation, such as C(r, t), 
C q (t), C(r, to), and C q {uj), and related versions containing 
more than one time or frequency arguments. Our con- 
vention for distinguishing them notationally, apart from 
displaying their arguments, is to use a tilde or bar to 
distinguish between the position and momentum repre- 
sentations, and a roman italic or calligraphic symbol to 
distinguish between the time and frequency representa- 
tions. 

A. Decay Function Fd(t) 

When the effect of interactions on the full Cooperon 
C(r = 0,t) is calculated within the influence functional 
approach, one naturally obtains results of the form 

C(0,t)~C°(0,t)e- F *W , F d (t) = ±<S e ff)rw . (2) 

Here C° is the bare Cooperon in the absence of inter- 
actions, and S e s is the so-called effective action. It is 
essentially the variance of the fluctuating difference of 
phases ^Sf and \Sb acquired while propagating along 

the two paths, S e s = 57- ({Sf — Sb) 2 )- In the case con- 
sidered here (linear coupling to Gaussian fluctuations), it 
is linear in the noise correlator (interaction propagator) 
and characterizes the effect of the environment on a pair 
of time- reversed trajectories whose interference gives rise 
to weak localization. Its average over all random walk 
trajectories yields the "decay function" Fd(t), describing 
the suppression of the Cooperon by decoherence. Hence 
the decoherence time t v = l/"f v can be defined by the 
conditional Fd{r v ) = 1. 



The decay function Fd{t) turns out [Sec. lIII D] to be of 
the following general form, for trajectories propagating 
during the time-interval [— i/2,t/2] [see Eq. ijSUjl]: 

I ft/2 rt/2 

F d (t) = Ta\ dt 3 dU 

11 J-t/2 J-t/2 
X J (dq) [P(q, |t 3 - U\) - P(q, |i 3 + U\)] 

xJme-W'-^iVV)* (3) 

It contains one time-integral for each of the two inter- 
fering trajectories. Besides, it is a product of a part 
describing the diffusive dynamics of the system under 
consideration (the second line) and the noise correlator 
of the effective environment (third line), which is inte- 
grated over all momentum and frequency transfers q,u>. 
In our notation, P(q, t') is the Fourier transform of the 
probability density P(r',t') for a random walk to cover 
the distance r' in the time t'. 

The fact that the second line of Eq. contains a 
difference between two rather similar expressions reflects 
the fact that the phases picked up along the two trajec- 
tories are related for fluctuations with sufficiently long 
wavelengths and/or low frequencies, and ensures that 
such fluctuations do not contribute to decoherence. The 
|*3 —ti\ and |*3 + *4| terms correspond to the "self-energy 
terms" and "vertex corrections" of diagrammatic calcul- 
tions of the Cooperon self-energy, in which the vertex 
terms are needed to cancel infrared divergences of fre- 
quency or momentum integrals. The simple and natu- 
ral way in which this cancellation arises in the influence 
functional approach is one of the latter's main advantages 
(the other being its physical transparency). 

To calculate P(r',t'), previous worka&S have usually 
averaged over unrestricted random walks (urw) that are 
not constrained to return to the origin, ignoring the 
fact that all paths contributing to weak localization are 
closed. We show in Sections IHI CI how P(t',t') and its 
Fourier transform P(q,t') may be calculated for closed 
random walks (crw) instead of unrestricted ones [Eqs. 
(|27J> and ail d in Section ITlI PI how the resulting 

more complicated decay functions may be evaluated. For 
d = 1 (but not for d = 2, 3), this improvement leads to a 
more accurate result for the numerical prefactor occuring 
in the decoherence time. [cf. Eq. (|44|) ]. The extent of 
the improvement obtained with the more accurate result, 
which is important for quantitative comparisons with ex- 
periment, is checked in Section IlII El by using it to cal- 
culate the magnctoconductivity for quasi 1-dimcnsional 
conductors with classical white Nyquist noise, and com- 
paring the result to the celebrated exact "Airy function 
expression" of AAKpS. [Using closed random walks also 
turns out to be an essential prerequisite for recovering 
the results of AAG from our theory, Sec. IVIES] . 

The difference between averaging over unrestricted ver- 
sus closed random walks can quite generally be summa- 
rized by the following formulas, found in Section IIII CI 
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(and confirmed in paper II): 



y/Dr v , given by 



Ficrw 
d 



771 urw 



(*) 
(*) 



C 1 ^ = 0,t) 



C°(r = 

r 1 



= 0,t) ' 

=o(*) • 



(4a) 
(4b) 



Here C 1 (r, t) is the first order term in an expansion of 
the full Cooperon C(r,t) in powers of the interaction, 
and Cg 7Jf _o(t) is its momentum Fourier transform in 
the absence of a magnetic field. We note that in both 
cases, Eq. {SJ represents the full Cooperon simply as a 
reexponentiated version of the first order term (either in 
momentum or real space) , but since the decay of the real- 
space Cooperon is required, the expansion is consistent 
to leading order in the interaction only if FJ rw (i) is used, 
which is why the latter gives more accurate results. 

In our paper, we successively present different versions 
of Eq. , which are distinct in the effective noise correla- 
tor (VV)° a of the environment (classical noise, quantum 
noise for single particle, or quantum noise for many-body 
situation with Pauli principle). They will all, however, be 
associated with some type of Nyquist noise, and factorize 
as Ti ( vv )°ql = ^0-' where Weff(w) will be called the 
"noise spectrum" . This will allow us [with some standard 
approximations, including an average over unrestricted 
random walks (urw)], to reduce Eq. (JSJ to the form 



,(*) ~ p d t 



dw 



rjl-d/2 



sin(wf) 



Uft 



(5) 



[the pd are given after Eq. 1|41|) ]. Note the presence of the 
infrared cutoff at u) ~ l/t that was mentioned above. 



9d(L v ) 



ha* 



(7i^r 1/2 , 71 = D^/h^y 



e*L%- d 



92, 



.92 = ha 2 /e 2 , 



(73TJ 1 / 2 , 73 = D(ha 3 /e^ 



(7) 



for d = 1,2,3, respectively. gd{L,p) conveniently lumps 
together all relevant material parameters into a single di- 
mcnsionless quantity, and will be used extensively below. 
Since good conductors are characterized by having a large 
dimcnsionlcss conductance, gd{L v ) 3> 1 [see Eq. <|18[) be- 
low], the last equality in Eq. 10 implies Tt v 3> 1. This 
means that for paths of duration t ~ r„, we have the 
inequality tT 1, which will be important below. 

As expected, we recover Eq. © when applying our 
influence functional approach to a single particle under 
the influence of classical white Nyquist noise in quasi- 1 
dimension: upon replacing W e ff(w) by W c i(a>) = T in 
Eq. (JSJ), we find 



f (t/r, 



AAKN3/2 

lp,l I 



(t/r v ,i) 



3/2 



(8) 



with c" rw = , thus the decay function is governed by 
the same decoherence time t** k as obtained by AAK. 
For the second equality, we used our convention of defin- 
ing the decoherence time via Fd(r Vy d) = 1, to absorb 
the numerical prcfactor into the decoherence time itself, 
yielding r V) i = T^ K (c" rw )~ 2/3 . If the calculation is done 
for closed random walks, the result is the same, except 
that the prefactor changes to c" w = \ph : /4. 



C. Quantum noise 



B. Classical white Nyquist noise 



If we consider a classical fluctuating potential V(x, t) 
acting on the electron (Sec. IHIjl . then (VV)^ is 



given by its symmetric noise correlator {VV)\ 



This 

was used in the seminal paper of Altshulcr, Aronov 
and KhmclnitskiiS, where they applied the classical 
fluctuation-dissipation theorem to obtain the thermal 
part of the Nyquist noise, for which W e ff(w) is simply 
given by the classical noise spectrum, W c i(w) = T, lead- 
ing to a decoherence rate that vanishes at T = (see also 
the semiclassical path integral analysis of Chakravarty 
and Schmidt and Stern, Aharonov and ImrjiiS). 

For example, for a quasi 1-dimcnsional disordered wire, 
AAK found^ 3 - 3 ! 



1 



k =7 v T = (TV7l 



2/3 



si(^T) 



(6) 



Here gd{L v ) is the dimcnsionlcss conductance for a quasi 
d-dimensional sample at the decoherence length L v — 



The case of a fully quantum-mechanical environment 
is more involved. If one considers the motion of a single 
electron in the presence of quantum noise (sqn) but 
in absence of a Fermi sea, one can apply the standard 
Feynman- Vernon influence functional approach^ (Sec. 
HV)l. In this way, one obtains Eq. © with {VV)\% 
replaced by the symmetrized quantum noise correlator 
(VV)f£ = \{{V,V}) qa - By detailed balance, this al- 
ways includes a factor coth(<D/2T) = 2n(Cu) + 1, (n be- 
ing the Bose function) and the resulting effective spec- 
trum turns out to be W S qn(w) = ^uj[2n(ui) + l], which 
describes both thermal (2n) and quantum (+1) fluctua- 
tions of the environment. (In this paper, temperature is 
measured in units of frequency, i.e. T stands for k^T/h 
throughout.) Physically, the quantum fluctuations incor- 
porate the decoherence due to spontaneous emission into 
the environment, which is possible even for an environ- 
ment at T = 0, if the single electron has a finite energy (in 
a metal, its energy is typically near ep). For such a single 
electron, quantum fluctuations thus lead to a finite T = 
decoherence rate (but for a physical model which is dis- 
tinct from the original one describing disordered metals, 
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which involve many electrons). The result for Fd tUIW (t) 
obtained for this type of noise turns out to coincides with 
the one obtained by Golubev and Zaikinifi, if, following 
them, we introduce an upper frequency cutoff by hand 
(to prevent an ultraviolet divergence), and take it to be 
the elastic scattering rate. 

However, diagrammatic calculations^^ [summarized 
in paper II, Section ??, see Eq. (??)] indicate that the 
presence of other electrons cannot be neglected, since the 
Pauli principle plays an important role in preserving the 
coherence of low-lying excitations in degenerate Fermi 
systems. Setting up a Dyson equation for the Cooperon 
in the momentum-frequency representation and extract- 
ing from the Cooperon self-energy the decoherence rate 
7^ for an electron with energy^ £ relative sf, one obtains 
a rate where the factor 2n(<D) + 1 is effectively replaced 
by 



2n{u>) + 1 + /(e + O) - f(s - Q) 



(9) 



In the literature, this factor often occurs in the form of 
the combination "coth(w/2T) + tanh[(e - w)/2T]". The 
Fermi functions ensure that processes which would vio- 
late the Pauli principle (u> 3> max[T, e]) do not occur, 
which turns out to eliminate the ultraviolet divergence 
mentioned above. However, in contrast to the influence 
functional approach it is rather difficult to properly in- 
clude vertex corrections in the diagrammatic approach. 
In fact, Fukuyama and Abrahams^ introduced a low- 
frequency cutoff 1/t^ by hand, which then has to be de- 
termined self-consistently. The neglect of vertex correc- 
tions also means that the decay function is always linear 
in t, whereas e.g. in quasi-lD the classical result is known 
to grow like t 3 / 2 (as emphasized by GZ in Ref.lTol). What 
is needed, evidently, is an expression for the decay func- 
tion that keeps both the vertex corrections and the Pauli 
principle (and thus is free from infrared and ultraviolet 
divergencies, respectively). This is the main goal of both 
papers I and II. 

In the present paper, we address the question of how 
to incorporate the Pauli principle in an influence func- 
tional approach. First we provide a heuristic discussion 
of decoherence in the presence of a Fermi sea (Sec. 0. If 
an initial perturbation creates a coherent superposition 
between two single-particle states A and A', the deco- 
herence rate (within a Golden Rule calculation, without 
vertex corrections) is given by a sum of a particle and a 
hole-scattering rate [see Eq. (|60[) ]: 

r„(A, a') = i[r e (A) + r fc (A') + r e (Y) + r h (\)} . (10) 

Inserting the usual scattering rates containing Fermi 
functions for Pauli blocking, one realizes that incorpo- 
rating the Pauli principle effectively means replacing the 
symmetrized quantum noise correlator by the following 
combination (see Eq. 1)64(1'): 



(W)l 



pp 



[f(e + Q)-f( e -B)]i(\y,V>]) 



(11) 



where we took the energies of the two relevant states to be 
nearly identical^, as they will be in a calculation of the 
zero-frequency conductivity. This formula corresponds 
to substituting for VV e ff(w) a "Pauli-principle-modified" 
spectrum W pp (uj), given by ^u) times Eq. and is con- 
sistent with the results obtained diagrammatically, e.g. 
by Fukuyama and Abrahams^. We then discuss the 
consequences of this comparatively simple prescription 
for adding the Pauli principle to an influence functional, 
and use it to calculate the energy dependence of j v . 



Sec. IVIII of the present paper I and all of paper II are 
devoted to a justification of this prescription from more 
rigorous approaches. In Sec. IV11I we demonstrate that 
our heuristic prescription yields a result that is equiva- 
lent to that recently obtained by one of us (JvD) by an 
analysis^ based on Golubev and Zaikin's exact influence 
functional for Fermi systems&2*i2*ii. Their expression for 
the effective action S e s contains Fermi functions that cor- 
rectly represent Pauli blocking. Indeed, it was pointed 
ouliSi that AAK's expressions for the decoherence rate 
can be recovered from this approach, by considering the 
action in momentum-frequency representation and prop- 
erly keeping recoil effects, i.e. the energy change e — > e^fu) 
that occurs each time the electron emits or absorbes a 
noise quantum. 

In paper II, we shall show how a diagrammatic anal- 
ysis based on an approximate solution of the full Bcthe- 
Salpeter equation for the Cooperon (including vertex cor- 
rections) leads to the same results as here. 



D. Results for decay function F^ p (t) 



The main novel results of our paper are contained in 
Sec. IVII where the decay function is evaluated explic- 
itly for the case of quantum Nyquist noise for an elec- 
tron moving in a Fermi sea of other electrons at thermal 
equilibrium. Using Eq. J2J) with the modified quantum 
noise correlator (VV)™ of Eq. JTTJ, wc find [Sec. NTK\ 
after averaging over the electron's energy that the de- 
cay functions F^ cxw (t) = (F^ rw {t)) E have the following 
forms (whose leading terms also follow from Eq. JSJ , with 
W pp (w) as spectrum) i2£ 



^f,crw W 
-^crwC*) 
^3,crw(*) 



(*/^,l) 3/2 



2 3/2 |C(i)l 1 



1 

\li(Tt) - (1 - TEuler) 



I - 



tt2 3 / 2 I 
3C(3/2) VTi 



, (12a) 
, (12b) 
(12c) 
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The leading terms depend on the "classical decoherence 



rates 



"34 



1 

T <pA 
1 

1 




3C(3/2) T 3 / 2 
r v ,3 " (7r325)V2 ^ 



(13a) 
(13b) 
(13c) 



which reproduce the results of AAK for classical white 
Nyquist noise (except that AAK's numerical prefac- 
tors are different, since our way of defining r^d is 
slightly different from theirs). The next-to- leading terms 
in Eqs. (|12fl generate the leading quantum corrections 
to these classical decoherence times. Extracting the 
modified decoherence times from the condition^ 

^Fcrw(w) = 1, we find 



T <P,2 = T <p,2 



1 + 



Iil(Tt Vi2 ) 



T v ,3 



T v ,3 



(14a) 
(14b) 
(14c) 



where h = 2 5 / 2 |C(±)|/(3^) = 0.8767, b 2 = 1 - TEuier = 
0.4228, and b 3 = tt2 3 / 2 /[3C(3/2)] = 1.134. Thus, the 
next-to-leading terms are parametrically smaller than the 
leading ones (confirming the conclusions of Ambegaokar 

and Vavilo\ii£) by g^ 1 ^ 2 (L v %) for d = 1, or l/ln,g2 
for d = 2, or g 3 1 (L v ^ 3 ) for d = 3. Our calculations 
therefore conclusively show that in the weak localization 
regime where gd(L Vi d) ^> 1, AAK's results for ob- 
tained by considering classical white Nyquist noise, re- 
main correct for quantum Nyquist noise acting on an 
electron moving inside a Fermi sea at thermal equilib- 
rium. Nevertheless, since it is not uncommon for weak 
localization experiments to reach the regime where the 
product Tt v is only on the order of 10 (e.g. Ref . l38l) . the 
corrections discussed here can amount to an appreciable 
effect [illustrated in Fig. below] . 

As a check of Eqs. (|12|l . we use them to calculate 
[Sec. IVI B) the first-order-in-interaction contribution to 

the weak localization magnetoconductivity, a^^ 1 ^, in 
the regime j v <C Jh <C T ', where 7# is the magnetic de- 
phasing rate. Reassuringly, this reproduces the leading 
and next-to-lcading terms of the corresponding results of 
AAGiii, obtained via an elaborate perturbative diagram- 
matic calculation, which keeps vertex corrections but is 
restricted to short times. We also show how to resolve 
an inconsistency between AAG's way of extracting the 



decoherence rate from a 



WL(1) 



and the results of AAK. 



Finally, we also discuss the energy-dependence of the 
decoherence rate f Sec. IVI fty . We calculate explicitly how 



the decoherence rate 7 yj d crosses over to essentially the 
energy relaxation rate cx e d ^ 2 as e is increased with re- 
spect to T, and find that the energy scale at which the 
crossover happens, namely Tg 2 ^ 3 \Lt), T\ng2 or T for 
d = 1,2 or 3, respectively, is parametrically larger than 
temperature for d = 1, 2. - This concludes our overview. 



III. COOPERON DECAY FOR CLASSICAL 
NOISE 

In this section we review how the decay of the 
Cooperon can be calculated using influence functionals 
for the case of classical noise. Although this is a stan- 
dard calculation, we shall cast it in a form that general- 
izes straightforwardly to the cases treated in subsequent 
sections, namely quantum noise [Sec. HVj and quantum 
noise plus Pauli principle [Sec. IV]. 



A. Definition of Cooperon 

The full Cooperon C(r,t) appearing in Eq. Q> for 
6aY^ can be written as a path integral 



C(r,t) 



r F {t/2)=r pr B (t/2)=0 

Vr F (t 3 ) / Vr B (t 3 ) 

r F (-t/2)=0 Jr B (-t/2)=r 

xA[r F (t a ),r B (t 3 )} (15) 



over pairs of electron paths with opposite start- and end- 
points, to be called forward and backward paths, with 
amplitude A[r F (t 3 ), r B (t 3 )]. The fact that they are time- 
reversed has been exploited to denote the start and end 
times of a path of duration t by ±t/2 (this yields time 
integrals over intervals symmetric around t = below, 
which turns out to be very convenient). Semiclassically, 
the path integral will be dominated by time-reversed 
pairs of diffusive paths, i.e. 



r F (t 3 ) = r(t 3 ) = r B (-t 3 ), 



(16) 



and for C(Q, t), these will have the same start- and end- 
points. 

In the absence of interactions and a magnetic 
field, the amplitude A[r F (■) , r B (■)} simply equals 
e *(s [r F (-)]-So[r B (-)])M wne re S [r(-)\ is the free action 
describing the propagation of a free electron through a 
disordered potential landscape. The corresponding free 
Cooperon propagator C° (r, t) is thus determined by the 
probability density for an unrestricted random walk (in 
c?-dimensions) to reach a volume element d 3 r separated 
from the initial point by a distance r, in time t: 

P urw (r,t) = -J-_ e - r2/4I)|t| (47rL>|i|)- d/2 . (17a) 

In the presence of a magnetic field H (which, for d = 1 
or 2, we shall assume to be perpendicular to the wire 
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or plane of the film), the free Cooperon is multiplied 
by a dephasing factor e~*/ Tff , where the magnetic de- 
phasing rate 7# = 1/th increases with increasing H 
(j H = 4DeH/(hc) for d = 3, or lH = D(eHa) 2 / (3c 2 H 2 ) 
for d = 1, 2, see Ref. 0). Thus, we have 



C°(r,t) = 0(t)P UIW (r,t)e~ t/TH 



(17b) 



[In contrast, the bare diffuson is magnetic-field indepen- 
dent: D°(r,t) = 6(t)P UIW (r,t).} 

Inserting Eqs. (0) and l|17b|) for the full Cooperon 
(7(0, t) into Eq. for the magnetoconductivity, the lat- 
ter can be written as 



r l+d/2 



dt e -*/Tff e --Fd(t) 

t 9d{L t ) 



(18) 



where L t — VDt. For H = 0, the integral is of or- 
der l/gd{L v ) (or larger for d = 2,3, since r e i/r v <C 1). 
(To see this, change variables to 2 = t/r^ and note that 
F(zr lp ) > 1 for z > 1.) Good conductors, which are char- 
acterized by the fact that the relative change 5aY h /ad 
in conductance due to weak localization is small even at 
zero magnetic field, therefore have l/gd{L^) <C 1. This 
is a well-known and very important small parameter in 
the theory of weak localization, which will be used re- 
peatedly below. (For d = 1, where it turns out that 
\/gd{L v ) — (71/T 1 ) 1 / 3 , this ceases to be a small param- 
eter at sufficiently small temperatures. This signals the 
onset of the regime of strong localization, which is beyond 
the scope of the present analysis.) 



B. Averaging over classical noise 

Let us now explore how the Cooperon is affected by 
interactions, or more generally, by noise fields. Generally 
speaking, these will cause the propagation amplitudes for 
the forward and backward paths to pick up random phase 
factors, hence destroying their constructive interference 
and causing the Cooperon to decay as function of time. 

In the case of purely classical noise, a single-particle 
description is exact, and the decay of the Cooperon can 
readily be evaluated using path integral a^ 2 ^ 5 . It is in- 
structive to review how this is done. Let us describe 
the noise, imagined to arise from some classical envi- 
ronmental bath, using a classical, real, scalar potential 
Vj = V(rj,tj), with correlator 



(dfc)e 8tec « {W)\ 



(19) 



(the superscript denotes classical, the prefactor — ih is 
conventional). Here we used the shorthand notation 
(dk) = (dq)(duj), with (dq) = d d q/[(2n) d a s - d ], (du)) = 
<&j/(27t), where k = (q,oj) is our standard notation to 
be used for momentum- and frequency-transfers between 



the electron and the bath, and kxij 
we abbreviate = r; — Tj , = ti 



- qrtj — Luttj , where 
tj (and, for future 



use, tij ~ ti+tj). The noise properties can be specified in 
terms of the Fourier components of the noise correlator, 
(VV)^. It is symmetric in q for homogeneous, isotropic 
samples. Moreover, for classical (but not quantum) noise, 
it is necessarily also symmetric in frequency, 



(20) 



because (ViVj)° is invariant under t^ — > tji. 

In the presence of a given configuration of the potential 
field Vj , the propagation amplitude for a pair of random 
forward and backward paths, r a (^3), with a = F/B, is 
multiplied by an extra phase factor ^(Sf-Sb)/^ ^ w ^ n 



ft/2 

(S F -S B ) = -i dt 3 V s a V 3a , (21) 

J-t/2 



a=F/B 



where Vj a = V(r a (tj), tj), and s a stands for Sp/g = ±1. 
The average of this phase factor over all configurations of 
the field Vj can be performed without any approximation 
if the field is assumed to have a Gaussian distribution^, 



,i{S F -S B )/n 



(22) 



where the "effective action" S c ff[r F (-),r B (■)] is a func- 
tional of the forward and backward paths and describes 
the effect of the environment on the propagating electron, 



1 f t/2 ^ 

5cff = 7TT7 / dt 3 dt± V S a S a > {VzaVia') , (23a) 
ltl J-t/2 



aa'=F/B 



(V 3a V 4a ,) = J (dk) e^te)— (t*)]-ot,«) (y a v a ,)_ . 

(23b) 

In the present section Mil (V a V a i)qa is (for all a, a') 
simply equal to the classical noise correlator (VV)^ of 
Eq. I|19fl . (The more general notation will become useful 
for reusing Eqs. (|23|l [and Eq. I|30|l below] in later sec- 
tions, which involve more complicated correlators.) Note 
that 5efi is purely real, because the classical correlators 
{VsaVia') are real. 

To obtain the effect of the environment on the 
Cooperon, e _Sotf /'' should be evaluated along and aver- 
aged over time-reversed pairs of paths [Eq. H16|) ]. The 
a = a' terms in Eq. (|23a|l then correspond to the "self- 
energy" contributions to the Cooperon decay rate, to 
adopt terminology that is commonly used in diagram- 
matic calculations of the Cooperon decay rate. These 
terms describe the decay of the individual propagation 
amplitudes for the forward or backward paths, corre- 
sponding to decay of the "retarded" or "advanced" prop- 
agators. The a a' terms in Eq. I|23a|l correspond to 
the "vertex corrections" to the Cooperon decay rate. 
The self-energy and vertex terms have opposite overall 
signs (sfSf = sbsb = 1 vs. sfsb = —1, respectively). 
Consequently, the contributions of fluctuations which are 
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Figure 1: Random walks from 1 = (ri,— 1/2) to 2 = 
(r2,+t/2), via 3 = (^3,43) and 4 = (r4,t4): For given times 
t, *3 and £4, the distance between and r± is overestimated 
when the closed random walk (left) is substituted by an unre- 
stricted random walk (right). This leads to increased decoher- 
ence and an underestimate of the weak-localization correction 
to the conductivity (see Fig. [5] dotted line). 

slower than the observation time (with uj < 1/t) and 
hence indistinguishable from a static random potential, 
mutually cancel, and do not contribute to decoherence. 
This is already apparent from Eq. I|21(l , even before aver- 
aging over Vf for sufficiently slow fluctuations, the two 
terms in i{Sp — Sb) cancel if r F (t 3 ) = r B (—t 3 ). 

C. Closed versus unrestricted random walks 

In order to explicitly evaluate the modification of the 
full Cooperon due to the fluctuating potential, we still 
have to average the factor e~ Soff / ri over diffusive paths 
r(t) (i.e. random walks) . An exact way of performing 
this average has been devised by AAK2^, by deriving and 
then solving a differential equation for the full Cooperon 
(which can be done exactly for the quasi 1-dimcnsional 
case in the presence of thermal white Nyquist noise). 
However, it is possible to obtain qualitatively equivalent 
results by a somewhat simpler approach (also used by 
GZ): Following Chakravarty and Schmid (CS)&, we ap- 
proximate the average over random walks by lifting the 
average into the exponent, 

(e- s ^ H ) ~e- F *« F d (t) = I (S cff ) rw . (24) 

\ / rw 

The "decay function" F(t) will turn out to grow with 
time (starting from F(0) = 0) and describes the decay 
of the Cooperon [cf. Eq. (??)]. By lifting the average 
into the exponent, we somewhat overestimate the decay 
of the Cooperon with time, since for any real variable x, 
the inequality (e~ x ) x > e~^ x holds independent of the 
distribution of x (Jensen's inequality) . 

At the corresponding point in their own work, CS 
make two further approximations when evaluating Fd(t): 
firstly, they do not evaluate the vertex corrections ex- 
plicitly, but instead mimic their effect by dropping (by 
hand) the contributions of frequency transfers Q < 1/t 
to the self-energy terms, i.e. they introduce a sharp in- 
frared cutoff in the latter's frequency integrals. Secondly, 



while averaging the correlators (V3 a V4 a ') cl of Eq. (I23bl) 
over random walks [i. e. averaging the Fourier exponents 
in Eq. (|23bf) ] . both CS and GZ employ the probability 
density P urw (r34, £34) for an unrestricted random walk 
to diffusively reach a volume clement d 3 r removed by a 
distance = r 3 — r±, in time ^34 1: 

/ e «?(r(t 3 )-r(t 4 ))\ = [ d 3 r 34 P urw {r 3i ,\t 34 \)e^ r3i 

\ I urw J 

= P urw (g,|t 34 |)-e- D<?2|tM| . (25) 

Here and below, position integrals like J dr 3 4 stand for 
a 3 ~ d J d d r; the prcfactor comes from the integral over 
the transverse directions, and it cancels the prefactor of 
P urw in Eq. (|T7aJ) . 

The two approximations discussed above are known 
to be adequate to correctly capture the functional de- 
pendence of the function F(t) on time, temperature, di- 
mcnsionless conductance, etc. In the following, however, 
we shall be more ambitious, and strive to evaluate the 
numerical prefactor of F(t) with reasonable accuracy, 
too. To this end, we have to go beyond the two ap- 
proximations of CS (dropping vertex terms and doing 
an unrestricted average), since both modify the numer- 
ical prefactor by a factor of order one: Firstly, we shall 
fully retain the vertex corrections; in effect, we thereby 
explicitly evaluate the actual shape of the effective in- 
frared cutoff function, instead of inserting a sharp cutoff 
by hand. Secondly, we shall perform the random walk av- 
erage more carefully than in Eq. 112511 . in that we consider 
only the actually relevant ensemble of closed*^ random 
walks of duration t, that are restricted to start and end 
at the same point in space: r(—t/2) = r(t/2) = 0. Thus, 
we use 

^ e if(r(ts)-r(t4)) ^ = J dr3 dr4 pcrw ^ ^ ^qr 3i 

= Pfo^feM). (26) 

Here P(™ t 12 )(3> 4) is the probability density for a closed 
random walk that starts at the space-time point (t^, £2) 
and ends at (r\,t\), to pass through two volume elements 
around the intermediate points (7-4, £4) and (r3,i3). For 
ti < ti < t 3 < t\ we have 

%™ tl2 )(3,4) (27) 
= f urw (ri3, *i 3 ) P mw (r 3 4, t 3i ) P urw (r 42 , t i2 ) 
P mw (ri2,t 12 ) 

The denominator ensures that the integral of 
pcrw over r4 ^ r3 yi e id s 1 [ as can be seen using 

/dr 3 P urw (ri3,ti 3 )P urw (r32,t32)_ = P urw (ri 2 , t 12 )\. 
We shall confirm below that p crw docs not depend 
on (ri,ti) and (7^, £2) separately, but only on the 
differences r\ 2 and t\ 2 , as anticipated on the left-hand 
side of Eq. l{27|. 

The probability density p crw obviously does not de- 
pend on the magnetic field. Note, though, that Eq. (|27|) 
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does not change if we multiply each factor p urw (r^ } t i} ■.) 
by a dephasing factor exp(— t^/ru) to obtain a bare 
Cooperon C° (fij ,Uj), since these factors completely can- 
cel out in Eq. (|T7J) . In the following, we shall thus use 
Eq. with P urw replaced by (7°, since this will be 
convenient when comparing to perturbativc expressions 
below, which are formulated in terms of C s. Performing 
the integrals J dr 3 dr 4 of Eq. ({26(1 by Fourier transforma- 
tion and using 



C° q (t)=9(t)e 



-E a t 



E„ 



Dq 2 + 1/th 



(28) 



for the momentum Fourier transform of the bare 
Cooperon C°(r, t) of Eq. H17b|) . we readily find: 



P, 



( C ™,ti 2 )(?'*34) 



J(dq) e^" C° q (t 13 ) C° q _ q (t 34 ) C° q (t 42 ) 
C°(r 12 ,t 12 ) 

e -D?t 34 (l-^) + iqr 12 ^ (29) 



In the limit t 12 — > oo this reduces to P ul ' w (q 7 t 34 ), as 
expected, provided t 34 is kept fixed. Note, though, that 
the latter condition is not appropriate for the evaluation 
of the long-time limit of the Cooperon, for which both 
time differences, ti 2 and £34, become large. 



D. General Form of the Decay Function F(t) 

Let us now evaluate the decay function F(t) = (S e ff) rw , 
starting from the effective action of Eq. ilL'Mli . Averaging 
the latter over time-reversed, random walks according to 
Eqs. l|16f) and (|25|) or l|2t)[) . the result can be written as 



F d (t) = - 



dU 



x(VV)f a 5P(q;T 34 J 34 ) 



dt 4 I (dq) / (duo) e 



(30) 



For classical noise, where (V a V a >}qo = (W)f Q , the 
"effective" environmental noise correlator appearing 
here likewise stands for the classical noise correlator, 
(VV)'~ = (VV)f Q . For the more general case that the 
(Va,Va>)qQ depend on a, a', as will be needed in our treat- 
ment of quantum noise below, the effective noise correla- 
tor is found to have the form 



(W)f Q 



(V F V F ) qa + (V B V B ) QQ 



<jLU 



(VbV f ) so + (V F V B ) 



(31a) 
(31b) 



where, looking ahead, we used the fact that the first and 
second lines are equal for all the types of noise to be 
considered in this paper. 

Eq. (|30|) also contains the object 



SP(q;T 34l T 34 ) = P(q,\t 34 \)-P(q,\t 3 



e -?Dtr 3i _ g - 



(32a) 
(32b) 



which describes the diffusive dynamics of the time- 
reversed trajectories. The first term in Eq. I)32|l arises 
from self-energy terms [with a = a' in Eq. (|23(l ]. the sec- 
ond from vertex corrections [with a ^ a']. Here T34 and 
T34 stand for 



T34 
T34 



= MA 

= |*34 1 A 



or 



or t 34 



T34 = [1 - \t 3i \/t] \t 34 \/t , (33a) 
[l-\t 34 \/t]\t 34 \/t, (33b) 



(£34 = t 3 + t 4 ) depending on whether the average over 
paths is performed over unrestricted or closed random 
walks [using P urw (g, |t 34 |) from Eq. (J2HJ) or P^ w t) (q,t 34 ) 
from Eq. H29(l ]. respectively. 

Since the time integrals in Eq. (|30|l are symmetric, only 



that part of (VV) e qQ that is symmetric under Co — > — Co 
contributes to Fd(t); moreover, if this symmetric part is 
real, so is Fd(t). 

The fact that the decay function is linear in interac- 
tion propagators has an important implication: when 
expanding both sides of Eq. (??) in powers of the in- 
teraction, the term linear in Fd(t) on the right-hand side 
of Eq. (??) must equal C 1 (0,t), the first-order contribu- 
tion to the full Cooperon (7(0, i), implying -i^ rw (t) = 
(^(O^VC^O,*) [cf. Eq. J£jJ)]. We added the superscript 
"crw" , because this relation turns out to hold only if the 
average over paths for Fd(t) is over closed random walks. 
The expression (??) for (7(0, t) thus amounts to a simple 
recxponcntiation of the first order interaction correction, 



(7(0, t) ~ C°(0,t) exp 



<5°(o,t) 



(34a) 



evaluated in the position-time representation (at r — 0). 
In contrast, if (following hitherto standard practice^) 
the average over paths is performed over wnrestricted 
random walks instead, i.e. if p urw is used instead of 
P crw in (HP, one finds -F$™(i) = J drC^ =0 (r,t) = 
Cl =lH=0 (t) [cf. Eq. l|iE|l]. In this case Eq. (??) yields 



C-(0,i)~(7 (0,i)exp[a^ =g=0 (t)] 



(34b) 



implying that here the first order correction in the mo- 
mentum-time representation is reexponentiated; conse- 
quently, the left- and right-hand sides of Eq. I|34b|) are 
not consistent when expanded to first order in the in- 
teraction. Hence, Eq. I|34a|) can be expected to be more 
accurate than (|34b|l . as will be confirmed below. 

To make further progress with the evaluation of F(t), 
we now exploit the fact that for all types of noise to 
be considered in this paper, the effective noise correlator 
factorizes into a frequency-dependent spectrum W c s(Co), 
symmetric in frequency, and a g-dependent denominator: 



vDq 2 



(35) 



This fact allows us to proceed quite far with the evalu- 
ation of Eq. (|30|) for Fd(t) without specifying the actual 
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form of W ff(d)) (which will be done in later sections): 
after changing integration variables to the sum and dif- 
ference times £34 and £34, Eq. 1)30(1 can be written as 

F d (t) = 2 / dt 34 W eS (t 34 ) / dt M 5P d (T M ,T M ) , (36) 
Jo Jo 



where the kernel 

W cS (t 3i )= I (duo) e-^W ctf (Q) 



(37) 



contains all information about the frequency-dependence 
of the noise correlator describes the range in time of the 
effective interaction, whereas the dimensionless quantity 



SPd(T 3 4,T 3 4) 



(dq) 



-q A DtT 3i 



-q'DtT 3i 



hvDq 2 



(38) 



is the same for all types of noise studied below, but de- 
pends on the dimensionality of the sample. Note that 
the integrand in Eq. (|38|) is well-behaved for small q 2 de- 
spite of the 1 /q 2 factor, because the self-energy and ver- 
tex contributions to SP(q; T34, T34) [Eq. i(32|) ] cancel each 
other for momentum transfers smaller than Dq 2 < 1/t, 
thus regularizing the divergence. This is precisely as ex- 
pected for density fluctuations with dispersion to ~ Dq 2 : 
on time scales of order t, fluctuations with frequencies 
below uj < 1/t appear to be essentially static, and hence 
do not contribute to decohcrence. 

The J(dq) and / dt 34 integrals in Eqs. ffljfy and ffify 
can be performed explicitly (see App. with a result 
for Fd(t) of the form 



Eq. I|39(l is a central result on which the following sec- 
tions rely. To find the decay function Fd(t) (and a corre- 
sponding decohcrence time), all that remains to be done 
is to determine the spectrum W c ff(k>) (which depends on 
the type of noise studied, and on whether the Pauli prin- 
ciple is taken care of), calculate its Fourier transform to 
obtain the kernel W c s(t 34 ) of Eq. 1(37(1 . and perform the 
integral / dt 34 in Eq. (|3T?|) . 

We close this section with a comment on the relation 
of the above approach to diagrammatic methods, e.g. 
the calculation of the Cooperon C 1 (0,i) to first order 
in the interaction (from which F c i(t) can be extracted). 
The key to our derivation of Eq. 139|) was to essentially 
work in the time domain, postponing any time inte- 
grals until after all momentum and frequency integrals 
had been performed, as exemplified by our definitions of 
P crw [Eq. ®, involving J(dq)}, SP d [Eq. involv- 
ing J{dq)] and W(t 34 ) [Eq. iflfo. involving f(du)]. In 
contrast, in diagrammatic calculations, the J dt 3 dt 4 inte- 
grals are performed first (namely when deriving the Feyn- 
man rules for the frequency- momentum), before any mo- 
mentum or frequency intetrals. However, for the present 
problem the resulting momentum integrals then take in- 
tractably complicated forms (see Add. !b|) . 

Nevertheless, for the case of unrestricted random 
walks, the leading asymptotic behavior of Fk )U rw(£) can 
be obtained rather simply by judiciously neglecting some 
terms that are sub-leading in the limit of large times. As 
shown in App. [5] this leads to the following expression, 



Fd,urvr(t) ^ Pdt 



W cS (to) 

duj 



Q2-4/2 



sin(wt) 



CO t 



(41) 



Fd{t) 



9d{L t ) J 



dt 34 W cS (t 34 )V d (t 34 /t) , (39) 



where Lt = V Dt. Explicit expressions for the func- 
tions Vd(z), which differ for closed or unrestricted ran- 
dom walks, are given in App. ^ However, it will turn 
out below that we really only need the leading terms in 
an expansion of Pd( z ) for small values of its arguments, 
because the decoherence rate is extracted from the long- 
time behavior of F d (t) ■ The leading and subleading terms 
of 7 7 " w for closed or p^ TW for unrestricted random walks 
(upper or lower entries, respectively) are given by: 



{i!£^)}-^ +0{z) > (40a) 



(,--) 



In 2: + 



2 + 0(z) 



(40b) 



r 3/2 











1 2 J 





1 + C(zlnz) 

+ 0(z 1 / 2 ). (40c) 



Thus, the difference between averaging over closed or un- 
restricted random walks turns out to matter for the lead- 
ing terms of only V\ (but not of ^2,3), implying', as we 
shall see below, that it matters for the leading long-time 
behavior of only F\(t) [but not of F 2j 3(t)]. 



with px = ^/27^/7^, p 2 = l/(g 2 2-K), p 3 = 1/( v / 2t 3 " it 2 ) 
[for 71,52,73, see Eq. Q]. This formula for F dlUTvr (t) 
is less accurate than Eq. I(39() . but perhaps physically 
more transparent, since formulated in the frequency do- 
main: The factor [1 — sin(u)i) / tot] acts as infrared cut- 
off, suppressing frequencies to < 1/t. Due to the factor 
Qd/2-2^ integral gets its dominant contribution from 
small frequencies of order 1/t for d = 1, gets logarithmic 
contributions from both small and large frequencies for 
d = 2, and is dominated by large frequencies for d = 3. 
In particular, for d — 2,3 an ultraviolet cutoff is needed 
at large frequencies to render the integral well-defined. 
As we shall see below, such a cutoff will be provided by 

Weff(w). 



E. Comparison with exact classical 1-D result 

To gauge quantitatively the difference between av- 
eraging over closed or unrestricted random walks, we 
shall now use the above results to calculate the mag- 
netoconductivity for a quasi 1-dimensional conductor 
and classical white Nyquist noise. For this particu- 
lar case, the average (e~ Self / n ) TTW was calculated ex- 
actly by AAK^ so that the resulting expression for the 
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Figure 2: The magnetoconductivity for a quasi 1-dimensional 
wire experiencing decoherence by classical white Nyquist 
noise, as function of t^ k /th, comparing AAK's exact Airy 
function result (14211 . solid line, against two approximations: 
the commonly employed average over unrestricted random 
walks in the path-integral exponent, [Eq. <34bH . and the 
improved version obtained by averaging over closed random 
walks [Eq. jgla) ]. resulting in Eq. (EU, with c? rw or cf w , re- 
spectively. For t v /th —* 0, the result for unrestricted (closed) 
random walks equals 0.808 (0.964) of the exact value. For 
t^ k /th 3> 1, all three curves become indistinguishable. 



magnetoconductivity 



:s:i 



WL 



nh [mA i(^-)]' 



(42) 



In 



can be used as benchmark for other approximations 
Eq. (|42|) . Ai is the Airy function, and r* AK is the deco- 
herence time of Eq. JSJ) . 

For the particular case of classical white Nyquist noise 
the noise correlator is given by T/vDq 2 , so that the ef- 
fective noise correlator of Eq. H30|) takes the form: 



i (W)£ - i {vv)% = 



T 



vDq 2 



(43) 



Thus, the weighting function in Eq. i|35|) is frequency- 
independent for classical white Nyquist noise, W c \(lo) = 
T, so that the corresponding Kernel is an infinitely sharp 
delta function, W c \(t^4) = TS^^). The J (H34 integral 
in Eq. (|39[l for Ff (t) is thus easily performed, yielding 



if(i) 



tT 



Cl 



(V^,i) 3/2 



1 

urw 



4 ' 
4 



(44) 



where 7% 1 = (ci^/txT) 2 / 3 . Depending on whether we 



average over closed or unrestricted random walks, two 



different values for the prefactor c\ are obtained. The 
decoherence time was obtained by solving F 1 (t v ,1) = 1 
for t Vi i, which reproduces Eq. ©, with an extra^I c\ in 
front of T. 

The fact that r£™ is somewhat larger than t£™ (the 
relative factor is (16/3tt) 2 / 3 = 1.423) can intuitively be 
understood from Fig.^ For a given time difference £34, 
the distance r 3 4 between the points r 3 and is overesti- 
mated on average when closed random walks are replaced 
by unrestricted random walks. Thus the latter give some- 
what too much weight to the effect of long wavelength or 
small momentum transfers, which dominate the integral 
in Eqs. (|38() . They hence somewhat overestimate the de- 
coherence rate, and consequently underestimate the de- 
coherence time and magnitude of the weak localization 
correction to the conductivity [cf. Fig. [2J. 

The weak localization contribution to the magnetocon- 
ductivity Sa^ 1 " of a quasi 1-dimensional wire can now be 
obtained by inserting Ff(t) of Eqs. 103} into Eq. (jT%)l . 
Fig. |3 compares the results so obtained using c" rw (dot- 
ted line) and c^ rw (dashed line) to AAK's Airy function 
result [Eq. H42f> . solid line]. Firstly and most importantly, 
all three approaches agree fully in their predition for r^,, 
which acts as the scale on which the magnetoconduc- 
tivity is suppressed as a function of increasing magnetic 
field (i.e. increasing 1/t#). However, the methods differ 
somewhat in their predictions for the magnetoconductiv- 
ity 8af h at B = 0, which gives the overall magnitude of 
the weak localization effect: Averaging over unrestricted 
random walks (dotted line) yields only qualitative agree- 
ment with the exact Airy-function result, deviating from 
it by about 20% at B = 0. In contrast, averaging over 
closed random walks (dashed line) gives excellent quan- 
titative agreement with the exact result, yielding prac- 
tically identical results for large and intermediate field 
strengths, with a maximal deviation of less than 4% in 
the limit of vanishing magnetic field. It is quite remark- 
able that such good agreement with an exact result can 
be obtained by means as elementary as the above. 



IV. COOPERON DECAY FOR QUANTUM 
NOISE WITHOUT PAULI PRINCIPLE 

Inspired by the simplicity and elegance of the above 
treatment of classical noise fields, we shall explore in this 
section to what extent a quantum bath can be similarly 
dealt with in a single-particle path integral picture: is 
it possible to identify suitably chosen "effective classical 
noise" correlators, such that the decay function F(t) is 
again of the form l|3l)|l . with (VV)f Q replaced by some 
suitably chosen effective noise correlator (VV)|g? The 
advantage of such a formulation would evidently be (i) 
that the results are sure to be free of infrared problems, 
and (ii) that the trajectory averages could be performed 
with the same ease as above. 

Of course, we know from the outset that a strategy 
based on mimicking quantum by classical noise fields 
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can never be exact or complete, because the correlator 
ZVV)q a of a quantum noise field differs from that of 
a classical noise field in an elementary but fundamen- 
tal way: In contrast to (V"V)SL , which is symmetric in 

frequency, (VV)^ is asymmetric in frequency, reflecting 
the asymmetry between energy absorption from (u> < 0) 
and emission into (w > 0) the bath, ui being the change 
in bath energy. In particular, the asymmetry manifests 
itself in the possibility of spontaneous emission events 
which are possible even at T = 0, and hence strongly 
affect the low-temperature behavior. Nevertheless, we 
shall see that when the effective action is evaluated along 
time-reversed paths as needed to describe the decay of the 
Cooperon, it is again governed by a noise correlator sym- 
metric in frequency, so that this decay can be described 
by a suitably-chosen "effective" classical noise field. 



A. Definition of quantum noise correlators 

We begin by recalling that for a free bosonic quantum 
field Vj = V(rj,tj), the two correlation functions 



ihl% = (ViVj) 



ihC? 



(VjVi) , (45) 



are not equal (as would be the case for a classical 
field), but are related, after Fourier transforming, by 
the detailed-balance relation Cf(ui) = e~" w Cq(u>). This 
implies that the symmetrized and anti-symmetric corre- 
lators, or, equivalently, the Keldysh, retarded and ad- 
vanced correlators 



-ihcZ = ({Vi,Vj}) 



R/A 



(46a) 

■iKC%" = ±9(±t ij )([V i ,V j }) , (46b) 
are related by the fluctuation-dissipation theorem, 



\iCf{u>) = lm£ R (u>) coth 



2T 



Im£«(|u>|)[2n(|u>|) + l] , 



(47a) 
(47b) 



where Im.£ R (cu) = -£fi{[V,V]) is an odd function of 
Cu. For example, the noise generated by the standard 
screened Coulomb interaction, namely quantum Nyquist 
noise (which we shall focus on for the remainder of this 
paper), can be described by taking 



(48) 



[see, e.g., Refs. fl5l orl2l|. Using the above relations, the 
quantum noise correlator can be written as 

\{VV). Q = hs£*(\Q\)[M\*\) + (49) 

In Eq. (|49[l . the contribution of n(|a>|) dominates at fre- 
quency transfers smaller than the temperature, for which 



the number of activated quanta is large (n > 1). The ad- 
ditional step function 6(u>), responsible for the asymme- 
try of (VV) , contributes even if the bath is at zero tem- 
perature (n = 0), and hence is sometimes said to reflect 
zero-point fluctuations of the bath. It describes the pos- 
sibility of spontaneous emission of energy by the electron 
into the bath, enabling excited electrons to relax to states 
of lower energy Of course, in a many-body situation, the 
rates for such relaxation processes will also contain Fermi 
functions that Pauli-block them if no empty final states 
are available. However, we shall defer a detailed discus- 
sion of Pauli blocking to SectionlVland completely ignore 
it in the present section, which thus applies only to sit- 
uations for which Pauli restrictions arc irrelevant. The 
latter would include a purely single-particle problem, or, 
for a many-body degenerate Fermi gas, an electron that is 
very highly excited above the Fermi surface, with plenty 
of empty states below to decay into. 



B. Averaging over quantum noise 

It is known^ that the effects of a quantum noise 
field V(r, t) on a quantum particle can be described 
in terms of classical (c-number) fields by proceeding as 
follows: one considers a path-integral with a Keldysh 
forward-backward contour, and includes the noise via 
phase factors e~ l $ dx i dt i v if/ n ( as part of e lSF / h ) and 
gijdxjdtjVjB/h ( as part f g-iSs/n) that contam two dif- 
ferent fields, VjF = Vp(rj,tj) and VjB = VB(fj,tj), 
on the forward- and backward-contour, respectively. [In 
the case of classical noise the Keldysh fields are equal, 
Vf = Vb = V.] The classical field correlators are related 
to the quantum noise field correlators (ViVj) of Eqs. (|45|) 
by time-ordering along the Keldysh contour: 



(ViFV JF ) 


ee (mv 3 ) = 


-W K 


+ C R 


+ £ A 


(v iB v jB ) 


EE (fViVj) = 




-C R 


-£ A 


(Vi B V jF ) 


= (Vi^) = 




+ C R 


-C A 


(ViFV jB ) 


= (V^i) = 




-C R 





C A h , (50c) 



The first set of relations in Eqs. I|5U|) follows from compar- 
ing the expansions generated when expanding the factors 
e — »/ dxjdtjVjF /h occur j n g j n a time-ordered path integral 
for the forward path, and e l f dx 3 dt 3 v jB /ft occuring in a 
backward path, with the corresponding expansions of the 
quantum time evolution operators Uf = Te~ l -> dx i dt 3 v i/ h 
and Ug = T e l I dx i dt i v i/ h occuring in Keldysh perturba- 
tion theory, respectively. The second set of relations in 
Eqs. (|50|l arc simply standard identities, following from 
the definitions © of Cf } ' R ' A . 

The effective action for a single-particle subject to 
such quantum noise (sqn), to be denoted by S^g", is 
obtained as for Eq. I|22|) : we now have to perform a 
Gaussian average over Vf , Vb , again with the action 
i(SF — Sb) of Eq. (|21fl in the exponent, but now with 
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Vj a = V a {r a (tj),tj), i.e. both the paths r a {t) and the 
fields V a differ on the forward and backward contours. 
The result for S^g" again has the form of Eqs. J22J, but 
now (V a V a '}- Q stands for the Fourier transform of the 

quantum correlators of Eqs. (|50() . In fact, e~ s <=ff l h is 
nothing but the Feynman- Vernon influence functional for 
a single particle interacting with a quantum bath. Note 
that in contrast to the effective action for classical noise, 
S^g 1 is in general complex, since the correlators (Vi a Vj a '} 
are all complex (because —i£f^j, £fj and Cfj are by con- 
struction purely real). 



C. Quantum noise spectrum without Pauli 
principle 

To describe the effect of quantum noise on the 
Cooperon, the influence functional e~ •>« ' has to be 
averaged over all time-reversed pairs of closed random 
walks. This can be done in the same way as as in Sec- 
tion EI3 The result for F s d qn (t) = ±(S c T>rw has the 
same form as Eq. <|3Ull . 



x(VV)f«5P(q;T 3i ,T 3i ), 



(51) 



where the effective noise correlator in Eq. ll.'jQl! now takes 
the form [obtained from Eqs. (|31J) and Eqs. I)5()|l]: 



(W)%~(W)'£ = -$i£K = 1({V,V}), Q - (52a) 



For the case of quantum Nyquist noise [Eq. 
jl(VV} s ^ can be written in the factorized form 
W S qn/ vDq 2 of Eqs. I|35|l . the corresponding spectrum be- 
ing 



W sqn (w) = i|w|[2n(|a)|) + 1] 



(52b) 



Just as for classical noise, this effective noise spectrum is 
symmetric in w and q and real [as follows from Eq. (|47b|) ]. 
which implies that F^ qn (t) is purely real. In other words, 
the imaginary part of the effective action S e g vanishes 
upon averaging over time-reversed paths; the reason is 
that (the Fourier transforms of) the purely imaginary 
T^i{£ R + £ A ) contributions from Eqs. (jijUaJ) and Ij50b|) 
cancel each other when inserted into Eq. I|31|) . as do the 
contributions ^\i(L R - C A ) from Eqs. l|5Uc|l and (|50d|) . 

The fact that (S^ n ) cr w is purely real along time- 
reversed paths has the following useful implication: for 
the particular purpose of calculating the Cooperon de- 
cay, it is possible to mimick the effect of a quantum- 
mechanical environment by a purely classical noise field, 
if we so wish, provided its noise correlator is postulated 
to be given precisely by (VV) S ^ of Eq. i.e. the 

symmetrized version of the asymmetric quantum noise 



correlator i|49[l . This can be verified by rewriting the ef- 
fective action in terms of even and odd combinations of 
V F / Bl namely 



v +i = h(y Fj + v Bj ) 



V-j = V Fj - V Bj . (53) 



It is then readily found that the decay of the Cooperon 
is governed only by the even field Vj + ; indeed, since their 
correlators are given by 



(V+iV +j ) (V+tV-j) 
(V-iV+i) (V-iV-j) 



= —ih 



lfK rR 
2 '"ij '"ij 

It o 



(54) 



SC(ii 



we see that the symmetrized noise correlator (VV) 
governing the Cooperon decay function F^ qn (t) is equal 
to the correlator (V+V+)qu of the even field V+j, whereas 
the correlators (V±V T )qu involving the odd field play no 
role in determining F^ qn (t). 



D. The effect of spontaneous emission 

It is instructive to analyse the differences between the 
classical spectrum W c i(w) = T and the quantum case 
Wsqn(w) of Eq. I|52b(l . Since both are symmetric in u>, 
the main qualitative difference between them is the pres- 
ence in the quantum case of spontaneous emission, lead- 
ing to an extra contribution that does not vanish at zero 
temperature. Although spontaneous emission only en- 
hances the scattering rate for transitions downward in 
energy, the preceding analysis shows that the asymmet- 
ric quantum noise spectrum may just as well be replaced 
by its symmetrized version. Physically, this is possible 
because both the upward and downward transitions are 
equally effective in contributing to dephasing (if they are 
allowed), and thus it is only their sum that matters for 
the dephasing rate. Schematically, we have 



(n + (n + l)) 
1, . 



2 ([n+- 2 ] + [n+-}), 



(55) 



where n = n(|u)|) is the Bose occupation number for the 
frequency transfer uj under consideration. 

This procedure is possible for a single, excited elec- 
tron without Fermi sea, or, in a many-body situation, 
for an electron so highly excited above the Fermi sea that 
Pauli restrictions on the available final states are negli- 
gible. In such a case, spontaneous emission processes 
evidently persist down to zero temperature and will thus 
cause the decoherence rate to remain finite even at T = 
[see also Ref.l2lj. Indeed, this can be seen explicitly from 
Eq. gTJl for -Fd q u n rw (i): replacing W e g(w) therein by the 
single-particle quantum noise spectrum at zero tempera- 
ture, wiqn -0 ^ (<I>) = \ \oj\, and introducing an upper cutoff 
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Jq U duj to regularize the ultraviolet divergence that then 

arises, one readily finds that ■F , JU w (i) = hpd&u^t, im- 
plying a finite decoherence rate at zero temperature: 



sqn(T=0) 



(56) 



V. DECOHERENCE AND THE PAULI 
PRINCIPLE 

The scenario discussed at the end of the previous sec- 
tion will of course change as soon as Pauli blocking be- 
comes relevant: consider a many-body situation, and a 
noise mode whose frequency Co is much larger than both 
the temperature and the excitation energy of the prop- 
agating electron. In such a case we expect that spon- 
taneous emission really would be severely inhibited by 
the lack of available final states. The present section 
is devoted to offering a heuristic understanding of these 
effects. (Previous, more formal, approaches'^ for deal- 
ing with Pauli blocking effects in a functional integral 
context are briefly reviewed in Section I VIII ) Remark- 
ably, we shall find that the decay of the Coopcron can 
again be described by a classical field with a symmetri- 
cal noise spectrum, but now containing an extra term to 
describe Pauli blocking, that turns out to block sponta- 
neous emission. Putting it differently, the Pauli blocking 
term counteracts the effects of the vacuum fluctuations 
of the environment. 



A. Early attempts to include the Pauli principle in 
path-integral calculations of decoherence 

The importance of Pauli blocking was recognized early 
on in the theory of decoherence in weak localization^. 
The simplest heuristic strategy to cope with this problem 
seems to be to derive the classical spectrum by apply- 
ing the classical fluctuation-dissipation theorem (FDT) 
to the linear response correlator. The 1 result is equiva- 
lent to replacing the [n(|u)|) + l] in Eq. fB^I by its low- 
frequency limit T/\lo\. This approach works well for the 
case of Johnson-Nyquist noised, which has a relatively 
large weight at low frequencies, so that these dominate 
anyway. Note, though, that more care has to be excer- 
ciscd for super-Ohmic baths, such as phonons. 

The question of how to include the Pauli principle has 
received surprisingly little attention in the early decoher- 
ence literature. The most concrete suggestion was due 
to Chakravarty and SchmidS^. For unexplained "gen- 
eral reasons" , i.e. probably in view of the perturbation- 
theoretic treatments^), they proposed the following re- 
placement as a way of incorporating the Pauli principle 
for the decoherence of thermally distributed electrons: 



2n(|w|) + l i ^ 2n(|t3|)+2/(|w|) 



(57) 



At low frequencies, \uj\ <C T, this yields the same factor 
T/\lo\ as the application of the classical FDT; moreover, 
it also provides an exponential cutoff at energy trans- 
fers ui larger than the temperature, thereby accounting 
for the absence of thermally excited bath modes at these 
high frequencies. The fact that Eq. I|57|) does not in- 
clude spontaneous emission at all (and actually vanishes 
at T — 0), is of no concern if we consider the deco- 
herence of an electron picked from a thermal distribu- 
tion, i.e. within T of the Fermi energy, for which spon- 
taneous emission would have been Pauli-blocked anyway. 
However, the absence of spontaneous emission would be 
a concern when describing highly-excited, non-thermal 
electrons, which, for a zero-temperature bath, would in- 
correctly be predicted not to relax at all. In other words, 
what is missing in Eqs. i|57|) is any reference to the energy 
e of the propagating electron. 

In the following sections we shall reanalyze these issues, 
but take care to include Pauli blocking throughout. Our 
conclusions turn out to qualitatively confirm the heuris- 
tic rule 157|) for thermal electrons, but are quantitatively 
more precise, and will also show how it should be gener- 
alized to deal with highly excited, non-thermal ones. 



B. Electron and hole decay rates 

In order to gain intuition about the effects of the Pauli 
principle, we shall first discuss the perturbativc calcula- 
tion of "golden rule" decoherence rates. Although this is 
only a first-order calculation in the interaction, the result 
is expected to be revealing nevertheless, since we know 
from Eq. (??) for F(t) that the decay function is needed 
only to linear order in the interaction, too. 

Consider a degenerate system of fcrmions under the in- 
fluence of a fluctuating environment that leads to transi- 
tions between the single-particle levels. The environment 
(e.g. a bath of harmonic oscillators) is described by a 
fluctuating potential V that couples to the fermions via 
some single-particle operator, and the correlator (VV)^ 
will determine the decoherence rate. For brevity of nota- 
tion, we do not consider a spatial dependence, V(r, t) — > 
V(t), and we shall assume the single-particle operator to 
connect any two levels with equal matrix element (which 
is set to unity in the following). The generalization to 
arbitrary coupling is straightforward. The golden rule 
decay rate for an electron in level A to be scattered into 
any other level is given by [Fig. |3j 

r> />oo 

r e (A) = — / (du>)(VV).[l-f(s x -u;)}D(e x -u>) 



(58) 



coth 
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i — ► coth 
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— tanh 
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2T 




2T 




2T 



sinh[M; 



(dCv) 2Im£ i? (w){ [l + n(u>)] [l - f(e x - w)] 
xD(e x -Cj)+ n(Q) [l - f(e x + u>)] D{e x + cj)} . 
where D(e) is the density of single-particle levels. In the 
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Figure 3: Sketch of two many-particle states forming a co- 
herent superposition, and two possible scattering processes 
contributing to T e (A) and I\(A), by emission and absorption, 
respectively. 



second equality, the first and second terms describe emis- 
sion and absorption processes, whereby the electron in 
the level A is scattered to a lower- or higher-lying empty 
level, respectively. At T = where n(u>) = 0, the only 
surviving process is spontaneous emission, and if e x ap- 
proaches the Fermi energy, e x — > 0, this process is sup- 
pressed, too, by the Fermi function [l — /(— w)] . 

Below we shall also need the rate for an initially empty 
state A to be filled, i.e. the decay rate of a hole, given by: 



r h (A) 



2?r 



(dQ)(VV) B f(e x + w)D(e x +w) (59) 

(du) 2Im£ i? (<I>)[ [1 + n{ui)] f(e x + Co) 
xD{e x + Co) + n{u>) f(e x - Q) D(e x - w)} . 



In the second equality, the first and second terms describe 
emission and absorption processes, whereby an electron 
from a higher- or lower-lying level is scattered into the 
empty hole level A, respectively. Again we see that at T = 
0, only spontaneous emission is possible, and if e x — > 0, 
this process is suppressed, too. 



C. Golden rule decay rates for coherent 
superpositions 

In order to calculate the decoherence rate, we have to 
consider a somewhat more complicated situation. Sup- 
pose we are interested in the linear response of the system 
to some perturbation (as is the case for the conductiv- 
ity calculation in the weak localization problem). If the 
perturbation scatters an electron from level A to level A', 
then the resulting state is of the form (1 + etpy'tp x )\^, 
where e <C 1 is small. Contributions of this type can 



occur if |<i>) has one electron in A but none in A'; 
apart from this restriction, |<I>) is some Slater deter- 
minant with arbitrary distribution of fermions over the 
other single-particle levels, and we will perform a ther- 
mal average over such states in the end. Effectively, 
we have thus created a coherent superposition of two 
many-particle states, which, for brevity, we shall call 
|1a,0a') + e|0>, 1a')- These are formed from an initial 
state with unoccupied A and A' by inserting a single ex- 



tra particle into a coherent superposition 



„t 



ec 



of 



these two levels (Fig. |3J- Our task is to calculate the 
decay rate of this coherent superposition, which under 
appropriate assumptions (discussed below) corresponds 
to the "decoherence rate" . Although we shall eventually 
need only the case e x = e x > (see Ref. |3^) , we shall for 
clarity distinguish the indices A and A' througout this 
subsection. 

The coherent superposition will be destroyed by any 
process that leads to a decay of one or the other many- 
particle state. This includes not only an electron leaving 
A or A', but also an electron entering the respective un- 
occupied state (A' or A). The total decay rate for the 
coherent superposition therefore is the sum of four con- 
tributions (after thermal averaging over the electron dis- 
tribution): 

v v {\ A') = I [r e (A) + v h {\') + r e (Y) + r A (A)] . (60) 

The first two terms give the decay rate for the state 
|1a,0a'), while the latter two refer to |0a, 1a')- The fac- 
tor | comes about because decoherence is due to the de- 
cay of wave functions rather than populations (the same 
is seen in usual master equation formulations of decoher- 
ence of systems with a discrete Hilbert space). 

In writing down Eq. (|6(J|) . we have assumed that all 
of the decay processes lead to decoherence. However, 
one may think of situations where Eq. H60(l would over- 
estimate the decoherence rate: For example, an electron 
traveling two different paths may scatter a phonon on 
both of these paths (similar to decay processes making 
the electron leave A and A'). But the interference is 
destroyed only if the wavelength of the phonon is suf- 
ficiently short to be able to distinguish the two paths 
from each other, since otherwise the information about 
the path of the electron is not revealed in the scattering 
process. In fact, disregarding this possibility amounts 
to neglecting vertex corrections in the diagrammatic cal- 
culation. Whether such an approximation is justified de- 
pends (among other things) on the operator whose expec- 
tation value is to be calculated in the end. This operator 
should connect the two states 1a, X >) and |0a,1a'), in 
order to be sensitive to the coherence of the state. It 
is necessary to specify this operator for each particular 
microscopic model, as well as the operator of the initial 
perturbation and the details of the system-bath coupling. 
However, such details are not important in the present 
section, since our aim here is merely to display the simple 
generic features of the golden rule decoherence rate in a 
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fcrmion system. 

It is illuminating to relate the calculation presented 
above to the decay of the single-particle retarded Green's 
function, G^(t), which appears in diagrammatic calcula- 
tions. According to the definition of Gjf, we have to 
consider the decay of the following overlaps: 



= 6{t) [(i>{u(t)^ \U(t)^ ) 

+(U(t)4> x ^o\4>xU(moj\ , (61) 



where |$o) denotes the ground state at T = (or some 
state over which a thermal average is to be performed). 
Here U(t) is the full time evolution operator, including 
the coupling to the environment. The first overlap can 
decay by two processes: Either the ket-state changes dur- 
ing the time t by the particle leaving the initial level A, 
or the bra-state changes by a particle entering A (before 
the time t). Thus, the decay rate for the first term is 



D. Pauli-blocked noise correlator (VV) PP 

\ / quj 



T e (X)+T h (X) 

cothf — ) + tanh 
\2TJ 



2tt 

ft 



(du>) Im£ K (w) x 

OO 



2T 



D(e x -u>). (62) 



By similar reasoning, the decay rate for the second term 
of Eq. I|61|) is given by the same expression. Note in par- 
ticular that a combination "coth + tanh" arises Eq. (|62|l ; 
this combination is well known from diagrammatic cal- 
culations of the decoherence rate in weak localization 
(see e.g. Ref. |29^> . At large positive energy transfers 
Cu 3> S\,T (emission of energy into the environment), 
this factor vanishes, due to the Pauli blocking of final 
states. Likewise, large negative energy transfers (absorp- 
tion of energy) are also forbidden, because the environ- 
ment does not contain thermal quanta needed for that 
process. 

Likewise, the decay rate of Gy{t) is found to be 
r e (A') + r^A'). Thus, (in the absence of vertex cor- 
rections) the decay rate for G^Gy coincides with the 
total decoherence rate T V (X, A') of Eq. (|6Ti|) . as expected. 
For the purpose of calculating the decoherence rate, for 
which we need the long-time behavior of the Coopcron 
(7(0, t) for Tt ^> 1, it suffices to take the electron and hole 
energies equals, E\ = ey , which will be done henceforth. 

In a related context, an explicit (non-diagrammatic) 
calculation of the decoherence rate has been performed 
for a fermionic ballistic Mach-Zehnder interferometer. 
This calculation is devoid of the complications introduced 
by impurity averaging but displays all the features re- 
garding the Pauli principle which have been discussed 
here^i. 



The discussion of decoherence and weak localization in 
Section ITTT1 had been purposefully restricted to situations 
where the Pauli principle played no role; in a many-body 
situation, this corresponds e.g. to a highly excited elec- 
tron state, with E\ so large that f{e\ ± Q) = 0. Taking 
the density of states D(e) to be constant henceforth, we 
then note from Eqs. f5%|) and (|5§|) that for / = the 
decoherence rate calculated in the previous subsection 
indeed depends on the symmetrized correlator only, 

r e (A) + r h (A) « [°° (dw) \{{V, V}) Q , (63) 

in agreement with the conclusion of Eqs. (|51|l and l|52|) . 

Moreover, we observe from Eqs. (|58|l and l|59[) that in 
the present context, "switching on" the Pauli principle 
(i.e. permitting / ^ 0) amounts to replacing \ {{V, V})_ 
in Eq. (jHSJl by 

V}) a + [/( £ a + - /( £ a - *)] !<[*>, V]) Q . (64) 

This yields a "Pauli-blocked" noise spectrum [compared 
to Eq. H52JI ]. which is, however, still symmetric and non- 
negative. 

Now, since the decay functions F(t) discussed in ear- 
lier sections and the golden rule analysis presented above 
are both linear in the interaction propagator, the lessons 
learnt from the golden rule about Pauli blocking should 
be directly relevant for F(t), too. Thus, we formulate the 
following hypothesis: if the first-order correction to the 
Cooperon C 1 were calculated by a proper many-body 
technique, the result for F(t) = -C l {Q,t)/C°{Q,t) will 
have the same form as Eq. I|51ll . 



dt 3 
, pp 



x(W)™ SP(q;r 34 ,T 34 ) 



dU / (dq) I {duj) e~ iat3i 

(65) 



except that the effective noise spectral function now has 
the following form, modified by the Pauli principle (pp): 



(vv) cS i — ► (vv) pp 



pp 



(jLO 



(66a) 



HiV, V}\ Q + [f(s + O) - f(e *)] \([V, V]) sa 



</UJ 



= Mm£f (o)){coth[^] + |(th_ - th + )}. (66b) 

Here e is the initial energy with which an elec- 
tron starts off on its diffusive trajectory, and th± = 
tanh[(e ± <D)/2T]. The combination "coth + tanh" aris- 
ing in Eq. I|66b|l is well known from diagrammatic calcu- 
lations of the decoherence rate in weak localization (see 
e.g. Ref.©. 

For the case of quantum Nyquist noise [Eq. (|48|) ]. 

can be written in the factorized form 
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Figure 4: The spectrum W pp (w) of Eq. 16711 as a function of 
frequency and electron energy. It has the properties W p p(0) = 
1, Wpp — > for Q > max{T,£}, and in the limit T 0, 
behaves as i|o)|5(|e| — \uj\). 



W pp {io)/vDq 2 of Eq. with a corresponding Pauli- 

principle-modified spectrum, 

Wpp(^) = ^{2n(M) + l + /( £ +M)-/(e-M)}, 

(67) 

shown in Fig. It has the very important prop- 

erty that it cuts off the contribution of all frequencies 
|c*> | > max{T, e}. In particular, for an electron at the 
Fermi surface (e = 0), the factor in curly brackets re- 
duces to the combination 2n(|u)|) + 2/(|a)|) anticipated 
by Chakravarty and Schmid [Eq. 1)57(1 ]. cutting off all 
frequencies \uj\ > T . At T = 0, the spectrum W pv (uj) re- 
duces to f |aS|0(|e| - \Q\). Moreover, at e = and T = 0, 
it yields {0 + 1 + 0—1} = 0, i.e. the new Pauli terms 
precisely cancel the spontaneous emission term discussed 
in Section TlV Dl as announced at the beginning of Sec- 
tion Thus, this spectrum can be expected to lead 
to a decoherence rate that vanishes for sufficiently small 
temperatures. We shall see below that this is indeed the 
case. 

The hypothesis that the use of (VV) 1 ^, in our formula 
for F(t) is the appropriate way to incorporate the Pauli 
principle into an influence functional approach will be 
shown to be correct in subsequent parts of this work: In 
Section IVII Al this replacement will be justified within 
the context of the functional integral analysis of decoher- 
ence of Ref. |3 Moreover, the Bethe-Salpeter analysis 
of paper II likewise turns out to lead to a decay func- 
tion [Eq. (II.??)] involving precisely the Pauli-principle- 
modificd weighting function of Eq. (|ulj)l : in particular, 
the diagrammatic calculation of C ll (0, t) performed there 
confirms explicitly that F(t) p d p = -C 1 (0,0/(7° (0, t). 



VI. RESULTS FOR DECAY FUNCTION F d pp (i) 

A. Energy-averaged Decay Function {F^ p (t)) e 

Since the correlator (VV) 1 ^ occuring in F pp (t) de- 
pends on the initial energy e with which an electron starts 
off on its diffusive trajectory, an average of the function 
e~ F d (*) over this energy still has to be performed, using 
the usual derivative of the Fermi function: 

(...)= = J de[-f(s)} ... . (68) 

We shall simplify the calculation by lifting the energy 
average into the exponent: 

(e-^Wj^e-^^, (69) 

thereby again somewhat overestimating the actual de- 
coherence rate [cf. discussion after Eq. I|24|) ]. The en- 
ergy average of the decay function, F^ p (t) = (F pp (t)) £ , 
has the same form as Eq. (I65|) . but now with an energy- 
averaged spectrum: 

W>)-<W)). = r(s!gy 2 . (70) 

which exponentially suppresses the contribution of all fre- 
quencies [ct> | > T, as anticipated above. The fact that 
Wpp (a)) < Wci(w) = T for all frequencies, but in particu- 
lar for u) > T, implies that the effective energy-averaged 
Pauli-blocked noise is somewhat less noisy than classical 
noise; thus, the decoherence times for the former can be 
expected to be somewhat longer than for the latter, as 
will indeed be found below. 

Evaluating the Fourier transform of Wpp (a)) (by clos- 
ing the J G?<I>-integral in Eq. I|37|) along a semicircular path 
in the complex plane), we find the kernel 

„9 , „ . . z coth z — 1 , 

W w (t 34 ) = irT 2 w(KTt 34 ), w(z) = — 2 , (71) 

sinn z 

where w(z) is a positive, peak-shaped function with 
weight J^° oo dzw(z) = 1. Thus, Wpp ^34) equals T times 
a broadened delta function of width l/T, and in the limit 
T — > 00, we recover Wpp(t 3i ) — > W c i(^34)- Inserting 
Eq. jnj into Eq. ffity. we obtain 




Since the decoherence time is defined by asking when 
F( T v) = an d 9d(Lip) 1 [cf. discussion after 

Eq. lfT5)l]. Eq. J72J) implies that Tt v > 1. Thus, to deter- 
mine F(t) for times t ~ r v , we may take the limit Tt 1 
in Eq. Q72[l . obtaining for the leading and next-to- leading 
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terms: 



ff (*) 
Ff P (t) 
Ff p (t) 



Tt 



9i(L t ) 
Tt 

52 (it) 
{Tt) 3 / 2 
33 (-^t) 



ci J= 

[02 ln(Ti) - 6a] 

b 3 " 



C3 



(73a) 
(73b) 
(73c) 



0.5826, 6§ rw 



-7Euler/(27T) = -0.09187, 6 C ,™ = ^Tr" 1 / 2 



0.06729, b%™ 



6 urw = 7r -3/2_ 



Here the prefactors c^ rw and c^ 1 ™ are the same as in 
Eq. 144(1 , reproducing the results obtained there for clas- 
sical white Nyquist noise. The prefactors c 2 = l/(2ir) 
and c 3 = J °° dzw(z)/(iry/z) = 3C(3/2)/\Ar325 = 0.2488 
are independent of whether the average over paths is per- 
formed over closed or unrestricted random walks, because 
the same is true for the leading terms of Vi and V3 in 
Eqs. I|4()|l. Hence, in contrast to d = 1, averaging over 
closed instead of unrestricted random walks yields no in- 
crease in accuracy for the leading terms of F pp (t) for 
d = 2 and 3. It does make a difference for the subleading 
terms, for which we find bf w = 6? rw = (27T)" 1 / 2 |C(|)| = 

1 - 7Eulcr)/(27T) 
crw 1 

3 - 2" ' > "3 

For unrestricted random walks, the leading terms of 
^durw(^) can a ^ so k e obtained with remarkable ease from 
its frequency representation (|41|) : replacing W c ff(w) by 
Wpp(<D) and evaluating the leading contributions to the 
integral in the limit Tt 1, one readily recovers the 
leading terms of Eqs. (|73|) [including the correct prefac- 
tors Cd\- 

Let us now calculate the full decoherence times r v d (in- 
cluding next-to-leading order corrections). For t ~ 
the next-to- leading terms in Eqs. (|73|l are parametri- 
cally parametrically smaller than the leading ones by 
g^ 1/2 (L vA ) for d = 1, or 1/ lng 2 for d = 2, or g^ 1/3 (L v _ 3 ) 
for d = 3. Therefore, we write r Vi d = T v .d(l + Sd), where 
5d <C 1 is a small correction induced by the next-to- 
leading terms, and first determine r Vi d, by setting bd — > 
and <$d — » 0. Then the condition^ F pp {r V) d) = 1 yields 
the the following selfconsistency relations and solutions, 



7v,i 

7^,2 
7^,3 



CiT 



3l(i V ,l) 

c 2 Tln (Tt Vi2 ) 
32(^,2) 
c 3 i 



7^,1 = (ciVtTT) 2 ^ (74a) 
c 2 T 

7^.2 = ln(5 2 /c 2 ) , (74b) 

32 



33 (L T ) 



7^,3 



c 3 T 3 / 2 



(74c) 



where 71, g 2 , 73 arc defined in Eq. 0, and Lt = y/D/T. 
These results reproduce those first derived by AAK for 
classical white Nyquis^jioise. They can be used to write 
the decay functions F pp (t) in the form (|12fl cited in the 
overview [Sec. Ill D] . 

Next we work out the corrections to the decoherence 
times due to the next-to- leading terms in Eqs. J73J). Re- 
instating bd ^ and solving the condition F pp (T v ,d) = 1 
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Figure 5: Comparison of the commonly employed leading ap- 
proximation for the decoherence time (thin lines) and the 
full decoherence time r v> d (thick lines), as a function of tem- 
perature, for dimensions d — 1 (full red lines), d — 2 (black, 
long dashed), d = 3 (blue, short dashed). The temperature 
T* has been (arbitrarily) chosen to give T*T ¥>) d(T*) = 10 in all 
cases (this amounts to fixing the values of the material param- 
eters 71, t/2 and 73). The magnitude of the correction is gov- 
erned by the small parameter [Tt^.J -1 , who's T-dependence 
is shown in the inset: In d = 1 (d = 3) the corrections become 
smaller (and the weak localization theory is applicable) to- 
wards higher (lower) temperatures, where [TY^e^T)] -1 <C 1. 
In d — 2, the correction amounts to a numerical constant fac- 
tor (shift on the log-scale, hardly visible in this figure). Note, 
in particular, the slow decay of the correction in the case of 



d = 1, where the correction falls off only like T 
Eqs. and <l?5aH . 



-1/6 



compare 



for &d{T v .d), we find 



Si = 



2 61/ci 
3# 



82 = 



6 2 /c 2 
ln(TT^ 2 ) 



(75a) 



leading to Eqs. (fH5|) for T Vt d- As anticipated above, the 
correction factors Sd are parametrically small, being of 

— 1/2 

order g 1 (£^,1) for d = 1, or l/lng 2 for d = 2, or 

— 1/3 

g 3 (Lif,3) for d = 3. We thus arrive at the most im- 
portant conclusion of this paper: the leading quantum 
corrections to the classical results for the decoherence 
rates and decay functions are parametrically small in the 
regime where weak localization theory is applicable. (The 
same qualitative conclusion was arrived at by Vavilov and 
Ambegaokar— several years ago by somewhat more indi- 
rect means). Nevertheless, note that the next-to- leading 
corrections are still parametrically larger than all futhcr 
subleading corrections, that could arise, e.g., from cal- 
culating F^ p (t) to second order in the interaction prop- 
agator, or from including crossterms between weak lo- 
calization and interaction corrections (as considered di- 
agrammatically in Ref. 0|)j smce such corrections are 
all smaller than the leading ones by at least l/gd(L v ). 
The leading and next-to-leading approximations to the 
decoherence times are plotted for all three dimensions 
in Fig. [S] As it is not uncommon for weak localization 
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experiments to reach to the regime where the product 



Tt v is only on the order of 10 (e.g. Ref. l38|) . we em- 
phasize that the corrections discussed here can amount 
to an appreciable effect. We also remark that, in d = 1, 
the relative size of the correction only falls off very slowly 
with increasing temperature (like T -1 / 6 ). 

The temperature dependence predicted for the cor- 
rected decoherence times t^^T) can be compared to 
experiment by proceeding as follows: Express F^t) of 
Eqs. I|73|> as a function of the parameters (i/r^d) and 
TT Vtd , by inserting T Vtd = r ip ^[l-(5 d (r ¥ , i( j)]. Calculate the 
magnctoconductivity cr^ L (H) numerically from Eq. (|T%)l . 
and for given T, adjust the parameter r^d such that the 
numerical curve as a function of magnetic field best fits 
the measured curve. Repeat for various T, and compare 
the function r ¥ , i d(T) obtained by this fitting procedure to 
the function r Vi d(l + <5d) predicted above. [If magnetic 
impurities are suspected to be present, insert a factor 
e -t/ T m m £ anc j treat the magnetic scattering 

time r m as a fit parameter. Spin-orbit scattering is not 
included in our analysis, but the corresponding general- 
ization should be straightforward.] 

To end this section, some remarks on the role of an 
ultraviolet cutoff seem to be in order at this point: for 
quantum noise in the absence of Pauli blocking, an ul- 
traviolet cutoff always has to be introduced to arrive at 
a finite result for the decoherence rate, to reqularize the 
contribution of spontaneous emission processes which oc- 
cur at all frequencies (see our discussion in Sec. II V D)| . 
In the full theory, Pauli blocking counteracts spontaneous 
emission and introduces via Wpp(w) an effective ultravi- 
olet cutoff at frequency transfers of order T. Remark- 
ably, for d = 1 (but not for d = 2, 3), the leading result 
for the decoherence rate can nevertheless be correctly 
obtained by simply employing the classical white noise 
spectrum W c i(<I>) = T (which contains no Pauli blocking, 
but no spontaneous emission either ) over the full fre- 
quency range up to arbitrary frequencies. The reason is 
that for d = 1, the dominant contribution to decoherence 
for time-reversed diffusive paths of duration t comes from 
frequencies Co ~ Dq 2 ~ 1/t [cf. Eq. (|38[) ]. which in the 
limit 1 <C Tt of present interest implies Co <C T; but for 
these frequencies, the spectrum yV pp (w) reduces simply 
to T, which equals the classical spectrum W c \(u>). 

In contrast, for d = 2,3 the large-frequency regime 
makes a logarithmic contribution for d — 2, and dom- 
inates for d = 3, requiring an ultraviolet cutoff to be 
present in the theory. AAK had to introduce such an 
ultraviolet cutoff by hand for the cases d = 2,3, be- 
cause they considered only classical white Nyquist noise, 
whose instantaneous kernel Wci(^34) = TS(t 3 4) involves 
no upper frequency cutoff, which is unphysical. Indeed, 
when one attempts to use it in Eq. (|39|l . the J dt 34 in- 
tegral would be ill-defined, since the 6^34) would pro- 
duce a ln(0) or 1/VO in Vi or V3 of Eqs. respec- 
tively. Equivalently, using W c \(Co) = T in Eq. (|4*T)) . 
the frequency integral would be ultraviolet-divergent for 
d = 2,3. AAK cured this problem by introducing, by 



hand, an ultraviolet cutoff at the scale of the temperature 
(taking \ q\ < y^T/D), and adding, by hand, a term to the 
decoherence rate describing the effect of electron-electron 
collisions with large energy transfers Co > T, which had 
been calculated earlier— within the framework of a ki- 
netic equation. (GZ's work implicitly questions this ap- 
proach, in that they used l/r e i (3> T) as upper cutoff.) 

Satisfactorily, an ultraviolet cutoff of precisely the type 
used by AAK arises automatically in our treatment of 
quantum Nyquist noise in combination with the Pauli 
principle, in the form of the energy-averaged spectrum 
Wpp(<D): it regularizes the large- frequency behavior of 
Eq. H41|) . without the need to consider processes with 
large energy transfers separately. Equivalently, in the 
time domain, it results in Wpp(i34) being a broadened 
(5(t34)-function of width 1/T. Thus, in a very natural 
(and perhaps somewhat more elegant) manner, we have 
confirmed the validity of AAK's use of temperature as 
an ultraviolet cutoff. Moreover, our explicit treatment 
of this cutoff was essential for accurately calculating the 
next-to-leading terms for the decay function and deco- 
herence times. 



B. Comparison with Magnetoconductivity of A AG 

It is instructive to check the use of the Pauli-blocked 
correlator (VV) PP introduced above against the results 

^ \ / qui ° 

of Ref. IT5I (AAG). There, the conductivity was calculated 
diagrammatically for the limit of a moderately strong 
magnetic field (for 7^ <C Jh "C T). In this regime the 
trajectories relevant for weak localization are so short 
(t < t h <C T v ) that the effects of interaction on weak 
localization are still small, so that it suffices to calcu- 
late 6a d A/h to first order in the interaction. At the same 
time, the condition 1/T <C th ensures that the premise 
for our calculations of t v in previous sections, namely 
1/T <C r„, still holds. AAG thus calculated the conduc- 
tivity diagrammatically to first order in the interaction 
and including all contributions of order 1/g 2 (including 
not only weak localization terms, but also interaction cor- 
rections and cross terms involving both). Among these 
1/g 2 terms, AAK identified the one that decreases most 
rapidly with magnetic field (largest power of th) as the 
one relevant for decoherence, and proposed to extract 
t v from it. They found that this term has the following 
form [Eq. (4.5) of Ref. l25l rewritten in terms of quantities 
introduced above], 



5a 



WL(1) 



-^[id* )m (d q) ^( { VV)Z) £ (76) 



C° q _ g (Q) \C° q (Co)\ 2 



[C°(0)] 2 C° q _- q {w) 



AAG evaluated the integrals using dimensional regular- 
ization, finding the following results^* [first two terms of 
their Eqs. (4.13>% we also evaluate their general result 
(4.11) for d — 3, sending (3 — d)^ 1 — > \u(th/t c \) in the 
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limit d — > 3, as appropriate for their dimensional rcgu- 
larization scheme! : 



r WL(1) 
1 , AAG 

r WL(1) 
0<Ji » ,V 
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^ttTth/2 
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3, AAG 
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3C(i 



(71-729)1/2 



lll(T_ff /t c i) 

8tt 3 VIVF 



03 



93(Lh)93{Lt) 



, (77a) 
(77b) 

(77c) 



Here . . . indicates sublcading terms with a weaker in- 
dependence. Moreover, AAG showed that the contribu- 
tion from cross terms between interaction corrections and 
weak localization (their Eq. (4.7) for ^0cwl) produce a 
Tff-dependence weaker than both the leading and next- 
to-leading terms of Eqs. fTTjl . 

Eq. 17tj|) can be reproduced from the analysis presented 
above, by calculating the magnetoconductivity using our 
first order expression for the Cooperon. To this end, 
use C 1 ' 6 = -C F% p civr [cf. Eq. g£J] in Eq. JTJ for the 
magnetoconductivity and average over e: 



end, a choice has to be made about the functional depen- 
dence on time of the full Cooperon C(0, t) (which AAG 
did not calculate explicitly), or, in our scheme, about the 
shape of the decay function F<i(t). Different choices for 
Fd(t) imply different "definitions" of r v , with different 
functional dependencies on temperature and magnetic 
field. In their Eq. (3.2), AAG chose to define l/r*^ G 
as a contribution to the "Cooperon mass" [in the sense 
of Eq. (II.??)], which implies that they assumed simple 
exponential decay for the Cooperon, e~^ 1 / Tff+1 / T ^ \ 
thus effectively making the choice F£ AG (t) = t/r AA d G . In- 
serting this into Eq. (|78a|l one finds 
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WL(l) 



r(2-|) t h 



n 1+d l 2 g d {L H ) 



(79) 



<p,d 



[reproducing^ A AG's (4.3)]. When equated to the lead- 
ing terms of Eqs. (|77(l . this yields for d = 1,2 



2gi(L H ) 



27T.92 



<p,i T > v,2 T\tl(Tth) 

reproducing^ AAG's (4.9), while for d = 3 we obtain 



(80a) 
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r l+d/2 



dtC°(0,t)(F™(t)). (78a) 
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V^25 g 3 (L T ) 
3C(| 
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(80b) 



dx 



x d / 2 g d (L H ) 



(78b) 



Now, for d — 3, Eq. (|80b|) reproduces the classical result 
of AAK [Eq. Ifflc))], However, for d=l,2 Eqs. pTajl for 
depend on magnetic field and hence are inconsis- 



(where L„ = y/Dr^). Now substitute Eq. ^ for tent with AAK's results [Eqs. CM] for r AA d K , which are 



^dcrw(*) m *° ^1- !|78a() . represent the <5P occuring 



therein via Eq. (|32a|l for closed random walks, with the P 
in Eq. (|32a|l standing for -P( C ™ (<7, £34), represented by the 
first line of Eqs. (|29() : then Fourier transform the three 
Coopcrons occuring in its numerator to the frequency do- 
main and finally perform all time integrals [Eq. (II.??) is 
helpful in this regard]; the result is found to be identical 

to Eq. H76f) . i- e - 60^^^ = ct^ag ■ The same conclu- 
sion can be reached by comparing AAG and our results 
after all necessary integrals have been performed: insert- 
ing Eqs. CZ2J) for FpO) into Eq. iffSb) for 5aJ L{1) and 
performing the time integral, we recover precisely AAG's 

results CHI> for SaJ^J . Thus, our theory is consistent 
with the calculation of AAG. Note, in particular, that 
the ncxt-to-leading terms of AAG's results for Sa^^j^ 
are also correctly reproduced in this manner; in our ap- 
proach, they are generated by the subleading_contribu- 
tions (the fr^-terms) to our decay functions F pp (t). This 
is a very useful consistency check. It illustrates firstly 
that our calculation of the next-to-leading corrections to 
the decoherence rate is correct, and secondly that the 
latter do not contain any contributions from the cross 
terms between weak localization and interaction correc- 
tions (which we did not calculate). 

AAG proposed to extract the decoherence times T Vt d 



magnetic- field independent, since AAK chose to define 
1 /t v as the decoherence rate which the Cooperon would 
have in the absence of a magnetic field. The reason for 
this inconsistency is that for d = 1,2, the Cooperon de- 
cay is not purely exponential in time [cf. Eqs. I|73|) ]. so 
that the usual strategy of adding inverse decay times to 
determine the total decay rate of two independent decay 
mechanisms cannot be used (as emphasized in the paper 
by AAK, after Eq. (32) of Ref. 

AAK's magnetic-field- independent results for the de- 
coherence time can be extracted from AAG's result for 
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tZ, AAG 



only if the correct functional form for the decay 



function is used. Indeed, inserting the leading terms of 
Eqs. for FfP(t) into Eq. fTHajl . we find 
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from their final results for 5a 



WL(1) 

d, AAG 



Eqs. If77fl]. To this 



When equated to Eqs. (|77jl . this yields 

t v ,i = (T^/4)- 1 , t v . 2 = 9 2/[Tc 2 ln( r *V>,a)], (82) 

implying 7^2 = ^7"l n (S2/c2) and thus reproducing 
AAK's results (with proper prcfactors included), as given 
by the rightmost equations of l|74|) . 
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To summarize this subsection, we conclude that, satis- 
fyingly, our results for the decay functions F pp (t) provide 
a bridge between the work of AAK and A AG: they allow 
AAK's results for r*^ K , obtained by treating classical 
white Nyquist noise nonperturbatively, thereby achiev- 
ing results free from infrared problems, to be extracted 
from AAG's results for Sa^^J , obtained by treating 
fully quantum noise perturbativcly, thereby incorporat- 
ing Pauli blocking and obtaining results free from ultravi- 
olet problems. The fact that our approach is able to make 
such a connection between two sets of established results, 
one nonperturbative but classical, the other quantum but 
perturbative, may be regarded as a strong indication that 
our method is fundamentally sound. 




C. Energy dependence of the decay function 

Instead of averaging over the energy he, it is also in- 
teresting to calculate the dependence of the decoherence 
rate 7^/ T on both temperature and energy, as would be 
relevant for an electron injected into a disordered metal 
with a definite energy (e.g. in a geometry such as that 
used by Pothier et a/4^). To the best of our knowledge, 
this energy-dependence has not been studied before. We 
shall now obtain it by analysing the energy-dependence 
of the decay function of Eq. (|65|> , for the case of closed 
random walks. 

To this end, we need the Fourier transform of the Pauli- 
principle-modified spectrum W pp (w) of Eq. I|67|l. which 
can be calculated by closing the / dui integral [Eq. 1)3 7|l] 
along a semicircular contour in the complex plane. The 
result can be written as W pv (t^4) = irT 2 w(et34 : , irTt^), 
where 



lip,! Hi 




Tfti 



cos y cosh z -\- (viz*) sin y sinh z — 1 

w(y,z) = y - y - . (83) 

2 sinh z 



Inserting this into Eq. I|39|) for Fit) we obtain the expres- 
sion 



Figure 6: (a) gi(L t )Fl P CIV (t) [Eq. JHD] as function of the 
parameters et and Tt. (b) The energy- and temperature- 
dependent decoherence rate 7>,i, defined from the condition 
•FFcrw( T ¥>,i) = 1, as function of T/ji and e/71. 



F "~ (i) =^P zm(w ^ w "(^)' (84) 



is given by the following expressions: 



which shows that F pp (t)gd(L t ) is a function of the pa- 
rameters irTt and x = e/T, or of irTt and et. Fig. HJa) 
shows the latter functional dependence for d = 1. 
Moreover, Fig. |(Jb) shows the corresponding energy- 
and temperature-dependent decoherence rate t Vi i, de- 
fined from the usual condition Ff p TW (T Vi i) = 1. 

To extract the decoherence rate analytically from 
F pp (t), we need its asymptotic behavior for large times. 
We find that the dominant behavior of Fd(t) for cither 
Tt ^> 1 or et ^S> 1 (or both), but arbitrary ratios of e/T, 



Tt 



9i(Lt) 

— I 

92{L t ) 

(Tt) 3 / 2 
93(Lt) 



Cl 



Fx{e/T) 



(Tt) 



1/2 



(85a) 



c 2 \n(Tt)+T 2 (e/T)}, (85b) 
1 



(85c) 



The crossover functions Td(ejT) are defined by the rela- 
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tions 



T d {x) = / dzw(zx/ir,z) l-[ln(z/ir) + 2]/tt } (86) 

i/(WI) 

for d = 1,2, 3, respectively. They have the properties 
f ^ d (0) + 0(a; 2 ) for s < 1 , 



<.//2 



for x > 1 



with ^i(O) = -0.6826, ^ 2 (0) = -0.1161, ^ 3 (0) = 0.2145 
and c 1 = a/2/tt, 5 2 = 1/(4tt), c 3 = 1/(3%/2tt 2 ). They gov- 
ern the crossover of the behavior of the decay functions 
from the regimes of small to large ratios e/T. 

For e/T -C 1, F d p {t) has precisely the same form as 
F pp (t) of Eqs. I)73|) . with the same prefactors ci and c 2 ; 
only the prefactor of F pp (t) is slightly different, namely 
c 3 = ^3(0) = 0.2145. Therefore, the decoherence rates 
have the same form l|74|l as that derived from F pp (t). 

In contrast, for e/T large enough that the large-x 
behavior of Td{x) dominates the behavior of F pp (t), 
Eqs. (|55|l reduce to 



(',1 



{et) d ' 2 



9d{L t ) 



t/ T V,d 1 lv,d 



9d(L E 



More explicitly, the corresponding decoherence rates and 
the crossover scales above which they apply are given by 

7^,1 = ci {lief 12 , for e > (T 4 /^) 1 / 3 = T 5l /3 (i T ) , 
7^,2 = c 2 e/g 2 , for e > Tin 52 , 

7 V ,3 = £ 3 (e 3 /73) 1/2 , fore>T. (89) 

Thus, for sufficiently high energies the decoherence rate 
has the same functional form as the inelastic energy re- 
laxation rate in quasi d dimensions 4 ^. It is interesting to 
note that for d = 1, 2 the crossover scales above which e 
has to lie for these relations to hold are parametrically 
much larger than temperature. 

The crossover behavior of r^d between the regimes 
of small and large e/T can be determined explicitly, 
if desired, from Eqs. I|85|) . using the usual relation 
F?(T v , d ) = 1. 

Note that for unrestricted random walks, the lead- 
ing terms of F% p (t) can again easily be obtained from 



will be quantitatively different from those of the more 
accurate Eq. (|84|l . However, qualitatively the crossover 
behavior is very similar. 

The fact that the functional dependence of r v i. 2 on 
e for e/T <C 1 is different than its dependence on T 
for e/T 3> 1, whereas for 7^3 it is the same, can be 
understood very nicely from the frequency representation 
(|4TT of -F 1 J > u rw : for T/e » 1, decoherence is dominated 
by high frequencies Q ~ T only for d — 3; for d = 2, 
the contribution from low frequencies of order Co ~ 1/t 
is as important as those from u) ~ T, and for d = 1, the 
contribution from Co ~ 1/t dominates. Thus, the infrared 
cutoff matters for d = 1 and 2, but not for d = 3. This 
is reflected in the fact that the first set of relations for 
7 Vj d in Eqs. lf74^l involve selfconsistency relations only 
for d = 1 and 2, but not for d = 3. In contrast, in the 
opposite regime of e/T 3> 1, decoherence is dominated by 
frequency transfers of order e not only for d = 3, but also 
for <i = 1, 2, so that the infrared cutoff is never important 
[to see this explicitly, use the T = version of W pp (u)), 
namely ^|cD|0(|e| — \u>\), in Eq. Q41|)]. Accordingly, none 
of the relations 7^ = Cde/gd(L e ) [Eq. JHSJl] involves a 
self-consistency condition, for every d. 



VII. RELATION TO THE WORK OF GOLUBEV 
AND ZAIKIN 



To close this paper, we shall now put the use of 
the Pauli-blocked noise correlator (VV) PP introduced in 
Eq. I|66|) on a firmer footing, bysummarizing how it fol- 
lows from the analysis of fiefs. and [3 111 R- e f- El Gol- 
ubev and Zaikin developed an influence functional for- 
malism and derived an effective action that explicitly 
and correctly included the Pauli principle. Indeed, a 
careful (if anfractuous) reanalysis of GZ's approach by 
von Delft24 has shown that one can, in fact, fully re- 
cover Keldysh diagrammatic perturbation theory from it 
(by starting from the initial, exact path integral expres- 
sion for the influence functional, and properly including 
fluctuations). However, when evaluating their effective 
action along time-reversed paths, GZ did not adequately 
account for the effects of recoil (as will be explained be- 
low). In Ref. I24I [see App. B.6.3], von Delft showed that 
the effects of recoil can be accommodated in the effective 
its frequency representation (fTTJl . by replacing W e ff(w) action by "dressing" the interaction correlators CfAuj) 



therein by W p 



(UJ 



[Eq. H67f) ]. which suppresses frequen- 



cies uj > max{T, e}. Evaluating the leading terms of 
Eq. (|4*T|) in the limits e/T — > and Tt^> 1, we recover 
the leading (e/T = 0) terms of Eqs. for F pp mw (t), 
with correct prefactors c d , and in the limits T/e — * 0, 
et S> 1, we recover Eqs. (JSSJ, with correct prefactors Cd- 
Since the derivation of Eq. I|41|) involved dropping some 
sublcading terms, the results which it will produce for 
the crossover behavior between the regimes of small and 
large ratios of e/T, and the subleading terms of F d (t), 



and Cfj (uj) (in the position-frequency representation) by 
suitably chosen "Pauli factors" 



[1 - 2/(e =F w)] = tanh[(e t w)/2T] = tff, 



(90) 



Instead of recapitulating the (lengthy) derivation of the 
latter conclusion, we shall now offer a plausibility argu- 
ment for it, based on the requirement of consistency with 
Keldysh perturbation theory 
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A. Describing Pauli blocking by dressed 
interaction propagators 

One way to see that the original Feynman- Vernon in- 
fluence functional e _Srff / fi . with S e s given by Eqs. I|23|) 
and (|50|l . cannot directly be used in a many-body situ- 
ation is that its expansion in powers of S c ff/h does not 
reproduce Keldysh perturbation theory, because the lat- 
ter contains Pauli factors, while the former does not. In 
Keldysh perturbation theory, the diagrams relevant for 
the calculation of the Cooperon have the property that 
each occurrence of an electron Keldysh Green's function 



th T [G R — G A ]ij(e =F ^) 



(91) 



B: 




j K 3 A 
(a) "type one" 



3 e+CO i 



R 00^ K 

A^4 




j K 

(b) "type two" 




with thzp = 1 — 2/(e ^w) = tanh[(e — Q)/2T], is always 
accompanied by either a retarded or an advanced (but 
never a Keldysh) interaction propagator attached to one 

of its two ends, C R J A (u)) or C R / A (uj). Since Pauli factors 

1 33 v ' it v ' 

enter in Keldysh perturbation theory only via G , this 
means that for the diagrams of present interest, every 
occurence of a Pauli factor is always accompanied by an 
C r / a (lu) propagator. To be consistent with this fact, 
the effective action in the influence functional approach 
must evidently contain the same combinations of Pauli 
factors and propagators C r / a (uj). This is not the case, 
however, for the C R l A {uj) occuring in Eqs. I|50|) . Thus, 
these propagagors have to be "dressed" by Pauli factors if 
we are to achieve consistency with Keldysh perturbation 
theory. 

The details of achieving consistency require two types 
of vertices to be distinguished (and keeping track of this 
strictly within the influence functional approach is the 
main technical accomplishment of Rcfs. .9:24): For ver- 
tices of "type one" [Fig. 0a)] , the arrows of the C R / A 
and G K correlators point in the same direction (i.e. both 
away from or both towards the same vertex), in which 
case Keldysh perturbation theory produces the combina- 
tions: 

C R 4p (w) Gf F4F (e - w) -> th_ C r 4f (w) , (92a) 
C A BAa , ip) G$ BiB (s-Q) ^ th_ C A Bial (Q) . (92b) 

For vertices of "type two" [Fig. 0(b)], the arrows point 
in opposite directions (one toward, the other away from 
the same vertex), which gives the combinations: 



C A Fi JCu)Gf F3F (e + Co) -> th+C 



3 F 4 r 



^ B ^)Gl jB (e + u>) 



(92c) 
(92d) 



The expressions on the right-hand sides of Eqs. (|92|l in- 
dicate how the L R ^ A t propagators in the effective action 

are to be dressed by Pauli factors (while £^4 , remains 
unmodified). Thus, the Pauli-modificd effective action 



Figure 7: (a) Vertices of "type one" and (b) of "type two" 
arising in Keldysh perturbation theory; the accompanying 
Keldysh Green's functions are G K (e^fu)), respectively, pro- 
ducing Pauli factors tanh[(e=F<Z>)/2T] that dress the associated 
interaction propagators C and C A [Eq. 1921 1. (c,d) Standard 
Keldysh rules generate diagrams such as both (c) and (d), but 
the latter vanishes upon impurity averaging, since it contains 
G A on the forward propagator. 



has the same form as Eq. (|23[1 . 

r t/2 

qpp _ J 



'1=4/ dhdh S * S a> (V 3 aV ia ,) PP , 

J-*/ 2 aa'=F/B 



(93a) 



(V 3a V 4a ,) pp = j (dk) e ^™~^ (V a V a ,) pp , 

(93b) 

and the correlators (V3 a V4 a ') pp have a form similar to 
Eqs. H5Ufl . but they are dressed according to the right- 
hand sides of Eqs. l|9"2|) : 



{VfV f ) pp = 






■f th_ 


effr) 


+ th + 




(VbV b ) pp = 


2 4 




- th + 




-th_ 




(V B V F ) PP = 


-¥ 




+■ th_ 




-th_ 




(V F V B ) PP = 


-¥ 




- th + 




+ th + 





(94) 

In this way, a single-particle influence functional formal- 
ism with suitably Pauli-dressed interaction propagators 
is able to mimick the essential features of the Keldysh 
many-body formalism. 

It should be emphasized that the possibility of us- 
ing such a dressing recipe is a feature peculiar to the 
Keldysh diagrammatics of disordered systems, as op- 
posed to a generic interacting many-fermion model (see 
Fig. 0: When inserting interaction lines into the two- 
particle propagator, standard Keldysh rules^i in general 
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also lead to contributions like G K {e — Cu)C K (to), which 
would spoil the above prescription of not dressing L . 
However, for disordered systems, such diagrams vanish 
to leading order, since instead of containing the usual re- 
tarded/advanced electron propagators G R I A on the for- 
ward/backward contour [as e.g. in Fig. Efc)], they con- 
tain G A / R there [as e.g. in Fig. |7Jc)], and after disorder 
averaging, 

The decay function F pp (t) = j-(S P g) rw , can be cal- 
culated in complete analogy to Sections IIV Bl and IIV CI 
Averaging as usual over time-reversed pairs of random 
walks using Eq. (|30fl . and calculating the effective noise 
correlator (VV)|? occuring therein using Eqs. I|31|) and 

Eqs. (|94|) . we now find that the dressed L R I A factors do 
not drop out (in contrast to what happens in Eq. i|52|) 
when calculating (VV) S ^ with undressed propagators). 
Instead, the effective noise correlator is now given by 
(W) e g% i-> {VV) v q l of Eq. JU, and the decay function 
(jH takes the form of F pp (t) of Eq. (EHJ. Thus, by using 
dressed propagators we recover precisely the Pauli-blocked 
correlator (VV)^ that we had conjectured on heuristic 
grounds in Section \V LA 

Further justification for using the latter noise correla- 
tor will be offered by the Bethe-Salpeter analysis of pa- 
per II, which leads to a decay function involving precisely 
the same Pauli-blocked noise correlator [cf. Eqs. (II.??)]. 
To be somewhat more precise: the Bethe-Salpeter anal- 
ysis of paper II improves upon the analysis of paper I 
by keeping track of the energy difference to (neglected in 
paper I, see Ref. 136^ between the energies e and e — lu 
of the particle and hole trajectories, which is why the 
effective interaction propagator Cf^^(oj) of Eq. (II.??) 
depends on lu. The latter propagator can actually be 
obtained from Eq. for (VV) c q % by the methods of 
the present section, namely by setting e — > (e — u>) in 
the Pauli-factors associated with the backward contour 
when calculating the dressed the correlators (V a V a ')qQ, 
i.e. in Eqs. i|92b[) and i|92djl . with corresponding changes 
in Eqs. ({94(1 . However, once the w-dependence is re- 
tained, the expression for Fd(t) needs another Fourier in- 
tegral J dcoe~ luJt , i.e. Eqs. (|30|l or (|65|l . and our intuitive 
determination of the factor 8P therein, become inade- 
quate. A more complete expression for Fd(t), including 
this additional integral, is given by Eq. (II.??). However, 
in the long-time limit max{T, e}t 3> 1 of present inter- 
est, for which lu/ max{T, e} <C 1, Eq. (II.??) reduces to 
Eq. (II.??), which reproduces our Eq. I|t)5|) . 

Finally, having been alerted to the necessity of keep- 
ing track of energy transfers, we note that the vertex dia- 
grams transfer energy between the forward and backward 
contours, thereby changing the energy of both. Strictly 
speaking, it is thus not fully correct to assign definite, 
fixed, energies e and e — ui to the upper and lower con- 
tours, since a succession of vertex insertions will produce 
an accumulation of frequency transfers. This effect is 
neglected in the influence functional approach. This is 
admissable, essentially because the dominant contribu- 



tion of vertex insertions comes from the infrared regime 
Cu ~ so that the (^-dependence in the remaining 

parts of the diagram may be dropped (they only gen- 
erate terms subleading in time). This point is discussed 
in more detail at the end of Sec. ?? of Paper II (and also 
in Sec. B.6.2 of Ref. Q). 



B. Importance of recoil 

To conclude this section, let us point out that the ef- 
fective action derived by GZ is essentially the same 
as our iSgg [Eq. i|93|l ]: the only difference is that in their 
approach, the dressed correlators emerge in the position- 
time representation of Eq. (|50|l [as opposed to the qtu 
representation of our Eqs. I|94|) ]. in which each C^ A is 
dressed by a Pauli factor [1 — 2/%, or [1 — 2(>\jj, where 
pij = {i(jj(t)ipi(t)) is the single-particle density matrix. 
They wrote their effective action as = ISr + Si, in 
which Sr contains Pauli factors and Si does not. When 
averaging over time-reversed pairs of trajectories, GZ 
used unrestricted random walks, F ( p z (t) = ■|(S'^ F z ) urw . 
Moreover, they employed a position- momentum repre- 
sentation in which they represented the Pauli factors as 
[l-2/(e(t))]. ^ 

Up to this point, their approach is essentially equiva- 
lent to ours. Differences arise at their next step: in order 
to perform the average over paths, GZ assumed the en- 
ergy of the diffusing electron to remain constant along its 
patb2£, arguing that its collisions with the static impu- 
rities are elastic, and hence replaced f(e(t)) by /(e). As 
shown by von Delft—, this assumption is equivalent to 
making the replacement e ± Cu — > e in the dressed prop- 
agators of our Eqs. 1)94(1 . although this is by no means 
obvious in GZ's formalism. In other words, they (unwit- 
tingly) neglected the ±lu in the Fermi functions occuring 
in the Pauli factors, and thereby neglected the recoil ex- 
perienced by the diffusing electron upon interacting with 
its environment and emitting or absorbing a noise quan- 
tum. [The fact that GZ neglect recoil was first pointed 
out by Erikson and HedegarcU&.] As a result, GZ's terms 
affected by Pauli blocking all cancel each other, mistak- 
enly causing (iS R ) urw to vanish, so that the resulting 
decay function F^ z (i) = ^(>S7) urw is identical to one 
for a single particle (no Pauli blocking) under the influ- 
ence of quantum noise. Indeed Fj" z {t) cani& be brought 
into the form of our ^"urwW [Eq. |5TJ]. As discussed in 
Section TlV Dl the resulting decoherence rate is the same 
as would be found for the physical situation of a single, 
highly excited electron moving through a disordered sam- 
ple very high above the Fermi surface, or in the absence 
of a Fermi sea, and interacting with a quantum environ- 
ment. Such an electron can lose its coherence even at 
zero temperature by spontaneous emission, uninhibited 
by Pauli blocking, which is why GZ obtained a nonvanish- 
ing decoherence rate at T = 0. Indeed, using il> u = 1/t c \ 
as upper cutoff in Eq. I|56|) for the decoherence rate of 
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a single particle experiencing quantum Nyquist noise at 
T = 0, we recover precisely GZ's zero-temperature deco- 
herence rate for all three values of d, including numerical 

r , sqn(T=0) GZ(T=0) r . „ rrn 

preiactors: 7„ d l v d \ l - e - our ^Q- <j5o|l re- 

produces the T = values of GZ's Eqs. (77) and (81) of 
Ref.U. 



VIII. CONCLUSION 

In this paper, we have shown why it is essential to 
incorporate the Pauli principle in a description of deco- 
herence by a quantum-mechanical environment, and how 
this can be achieved in an influence functional description 
of weak localization. We have explicitly demonstrated 
how Pauli blocking counteracts the effects of spontaneous 
emission to ensure that the decoherence rate vanishes for 
sufficiently small temperatures and energies. 

At the beginning of this paper, we offered a review 
of the influence functional approach for a single electron 
in a classical environment, pointing out along the way 
that quantitative improvements can be achieved by per- 
forming trajectory averages with respect to closed (as 
opposed to unrestricted) random walks SL. We then ex- 
plained how to extend the influence functional method- 
ology to the case of a quantum environment, taking due 
account of the Pauli principle. 

In our approach, a fully quantum-mechanical environ- 
ment may be treated in complete analogy to the much 
simpler case of classical noise. To this end, the quan- 
tum noise spectrum is modified in a well-defined way 
that involves Fermi functions, in order to take the place 
of the classical noise spectrum in the calculation of the 
Coopcron decay function (and the resulting decoherence 
time). We have shown how this replacement can be mo- 
tivated heuristically using a transparent physical picture 
for the decay of a superposition of two many-body states, 
which shows that the basic idea of the present approach 
is general enough to be extended to situations different 
from weak localization. In limits where the Pauli prin- 
ciple is ineffective, our theory reduces to the Feynman- 
Vernon influence functional, with a nonvanishing deco- 
herence rate even at zero temperature. In contrast, for 
electrons propagating near the Fermi surface Pauli block- 
ing is very important: it serves to essentially suppress 
the decohering effects of quantum fluctuations, confirm- 
ing the results of Altshulcr, Aaronov and Khmclnitskii, 
which were derived by keeping only the classical (ther- 
mal) part of the fluctuations. Moreover, our approach 
has also enabled us to calculate quantitatively the leading 
corrections to the decoherence rate, and to discuss in de- 
tail the energy-dependence of the decoherence rate (and 
the Cooperon decay function), also for energies higher 
than the temperature. 

Paper II will be devoted to substantiating these re- 
sults with the full machinery of Kcldysh diagrammatic 
perturbation theory. 
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Appendix A: EVALUATION OF P(z) 

In this appendix we derive Eqs. I|39fl for Fd(t) and (|40|) 
for Vd(t)i starting from Eqs. (|36() . 

The latter contains the function <5Pd(734, "?34), defined 
in Eq. J38J)), which may be calculated as follows: We 
render the momentum integral in Eq. I|38[) Gaussian using 
the integral representation 



q 2 D 



q 2 D 



tr I dxe~? DtTX , (Al) 

1/t 



then perform the Gaussian J (dq) integral, and finally the 
auxiliary J dx integrals, obtaining (with L t = \J Dt): 



-,2-d 



SPd(T, T) 



7T d / 2 9d(L t ) 



\ ln(?/r) \ for \ d = 2 



_L 1_ 



d = 3 
(A2) 



Inserting these expressions into Eq. (|36[) and performing 
the 34 integral results in Eq. 1|39|) for the decay function 
Fd(t), where for closed random walks we obtain 



r l/2 



{2z-S)[(l-z)z\ 2 + icos" 1 (2z-l) 



1 



7> 2 crw (z) = - 2(2-1) -In 2 



(A3a) 



1 



r 3/2 



[{l-z)/z} 2 - cos" 1 (22-1) 



(with cos 1 (2z — 1) g [tt,0] for z G [0,1]), whereas for 



2G 



unrestricted random walks we get 
1 



?>r w (z) = 

(z) = 



■purw 



7rV2 
1 
TT 



p™(z) 



r 3/2 



f(l-«)* -4(l-z)z5 
(l-z)(ln[(l-z)/z] -l) , (A3b) 
(1 -;»)/«* -2(1- *)l 



Expanding Pd{z) for small values of its argument yields 
Eqs. lj4*0)l of the main text. 



Appendix B: FREQUENCY REPRESENTATION 

FOR F d (t) 

In this appendix we discuss the derivation of Eq. 1)410. 
which expresses -Fk,urw(i) in terms of a frequency integral. 
We shall begin, however, more generally, by considering 
the case of closed random walks, starting from Eq. i|3U|) 
for F diCrw (t), with (VV}~ j given by Eq. §J$. Setting 
i 4 — > — i 4 in the vertex term of Eq. I|32a[) for SP CTW yields 
the combination [ e - ?lSt34 - e~^™]P?™Jq, i 34 ). Express- 
ing the latter factor through the first line of Eqs. (|29[1 . 
we obtain 



(dq)(dq)(du> 



Weff(^) 



(Bl) 



TwDq 2 C°(0,t) 
where we have defined the kernel [with A = (E q _ q — E q )\ 



K-- = 2 

qqu 



2e 



-E a t 



«3 

dti 

2 

At 



-E a t~At 3 



sin(wi) 



(B2) 



+ 



e ~t(\+i.a) _ 1 
(A + wl>) 2 



-tA 



The corresponding expressions for F^™(t) have the same 
form, but [in accord with Eqs. (@}] without the / (dq) 
integral and the factor 1/C°, and with E q = in the 
integrand, so that A ^ A = Dq 2 , K qqa /C a ^ K^. 

To proceed further, the spectrum VV e ff(w) has to be 
specified. For the case of classical white Nyquist noise, 



where it is given by W c i(w) = T, the frequency integrals 
can be performed explicitly by contour methods, yielding 
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For F^ urw (t), this result stems purely from the first line 
of Eq. (|B2|I . since the two terms in its second line each 
(separately) give zero. [Note that Eqs. (|B3|I are actu- 
ally obtained most easily by performing the J(du)) inte- 
gral in Eq. (|B1|I before the time integrals.] Both expres- 
sions for F£ l (t) are free of infrared divergences as q — > 
[in accord with the discussion after Eq. 1)38(1 ]. but for 
d = 2,3, they are ultraviolet divergent, as discussed at 
the end of Sec. I VI Al For d = 1, the momentum integrals 
can be performed by first rendering them Gaussian, us- 
ing the integral representation (1 — e~ b )/b = J Q dxe~ bx , 
and performing the auxiliary dx integrals last, whereupon 
Eqs. (|44J) are recovered. 

For non-trivial choices of W c ff(w), such as W sqn or 
Wpp for quantum noise, the J (dw) integral has to be 
performed last. To make progress with the momentum 
integrals, let us make some simplifications (which the 
main text manages to avoid): Firstly, we shall hence- 
forth study only the case of unrestricted random walks, 
Fd,urw(t) [whose asymptotic large-i behavior differs from 
that of Fd^ C rw(t) only for d = 1, as shown in the main 
text] . Secondly, we shall retain only the term in the first 



line of Eq. l|B2f) for K\ qQ , which dominates the behavior 
for large t, and drop the "sublcading" terms in the sec- 
ond line (which yield zero for W c i, as mentioned above). 
Performing the J (dq) integral in Eq. (|B1J) using 



(dq) 
tt J hvDq 2 



Ac 



d/2-2 



(B4) 



we obtain Eq. I|41|l of the main text, with pi = s/Hji/ir, 
P2 = 1/(32271"), p 3 = l/iV^fyiT 2 ) [for 71,32,73, see 
Eq. 0]. 
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